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ABSTRACT 

This  work  was  undertaken  to  determine  the  pressure  drops  in  the 
oxidizer  flow  passages  of  a  compact  dual -combust ion  liquid  rocket 
engine  by  means  of  fluid  dynamic  model  tests,  using  air  instead  of 
liquids,  to  establish  design  data.   Investigated  also  were  flow  dis- 
tributions and  flow  instabilities  by  measurements  and  flow  visuali- 
zation methods.   Changes  are  indicated  to  improve  the  original  design. 

The  air  tests  were  carried  out  over  a  range  of  Reynolds  numbers 
to  establish  the  effect  of  compressibility  on  the  pressure  drop  loss 
coefficients.  Water  tests  were  undertaken  to  check  the  validity  of 
the  results  and  of  the  proposed  calculating  method. 

The   tests  were  conducted  at   the  Turbo-Propulsion  Laboratory  of 
the  Naval   Postgraduate   School,   Monterey,   California. 


Thesis  by  Robert  Joseph  Kelly  entitled:  "An  investigation 
of  the  Oxidizer  Flow  Passages  in  a  Dual-Combustion  Liquid 
Rocket  Engine" 
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SECTION  1 
INTRODUCTION 

The  science  of  astronautics  has  progressed  at  an  exponential  rate 
in  the  last  decade.   The  development  of  liquid  rocket  engines  capable 
of  producing  large  thrusts  at  high  specific  impulse  has  been  responsi- 
ble for  this  progress.   These  engines  require  very  large  flow  rates  of 
fuel  and  oxidizer.   To  keep  their  size  and  weight  at  a  minimum  the  com- 
bustion pressures  are  high,  limited  only  by  structural  and  turbo-pump 
design  considerations.   The  flow  passages  inside  the  units  are  intri- 
cate and  costly  to  manufacture.   Pressure  losses  in  these  passages  are 
high  and  must  be  known  in  order  to  establish  pump  and  turbine  design 
point  data.   Accurate  prediction  of  these  losses  by  analytical  means  is 
not  possible  except  for  flow  passages  with  simple  geometries.   There- 
fore, it  is  advisable  to  conduct  model  tests  with  fluids  at  low  pres- 
sure, preferably  air,  to  determine  pressure  loss  coefficients  experi- 
mentally which  can  then  be  used  to  predict  the  pressure  losses  in  the 
actual  engine. 

The  ARES  (Advanced  Rocket  Engine  Storable)  propulsion  unit  is  a 
compact  dual -combust ion  liquid  rocket  engine  designed  by  the  Aerojet- 
General  Corporation,  Sacramento,  California.   Professor  M.  H.  Vavra  of 
the  Department  of  Aeronautics,  Naval  Postgraduate  School,  Monterey, 
California,  has  designed  full-scale  fluid  dynamics  models  of  this  unit 
and  has  conducted  preliminary  tests  of  the  oxidizer  flow  passages  with 
air.   The  models  have  been  made  available  to  the  writer  by  Aerojet- 
General  and  were  used  to  carry  out  more  extensive  and  more  detailed 
testing  with  improved  instrumentation. 
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The  purpose  of  this  thesis  is  threefold.   First,  the  air  test 
results  are  used  to  predict  the  pressure  drops  in  the  flow  passages, 
flow  distribution  within  the  unit,  and  velocity  distribution  at  the 
entrance  to  the  combustion  chamber  for  a  liquid  flow.   Second,  the 
effects  of  compressibility  and  Reynolds  number  on  the  test  results  are 
investigated,  and  the  validity  of  the  testing  method  is  established. 
Third,  changes  in  design  are  recommended  based  on  the  test  results  and 
flow  visualizations. 

The  author  wishes  to  acknowledge  the  guidance  provided  by  Profes- 
sor Vavra  and  extends  sincere  thanks  to  Mr.  R.  W.  Savage  for  his  assis- 
tance in  performing  experiments.   Acknowledgment  is  due  also  to  Mr.  R. 
Beichel,  Aerojet-General  Corporation,  for  the  loan  of  the  models  and 
his  endeavors  to  apply  model  testing  techniques  to  the  design  of  ad- 
vanced rocket  engines. 
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SECTION  2 
DESCRIPTION  OF  ARES  UNIT  AND  TEST  MODEL 

A  three-dimensional  schematic  representation  of  the  oxidizer 
passages  of  the  ARES  propulsion  unit  is  shown  in  Fig.  1.   The  oxidizer 
is  composed  of  a  mixture  of  nitrogen  tetroxide  and  hydrazine  and  it  is 
used  for  regenerative  cooling  of  the  combustion  chamber  and  the  rocket 
nozzle  before  being  mixed  with  the  fuel.   Oxidizer  enters  the  outer 
annulus  through  a  centrifugal  pump.   This  annulus  contains  turning 
vanes,  not  shown  in  Fig.  1,  to  distribute  the  flow  uniformly  to  the 
tubes  that  form  the  outer  skin  of  the  nozzle.   Each  tube  is  an  inte- 
gral unit  consisting  of  a  down-flow  and  an  up- flow  passage.   The  en- 
trance and  exit  holes  are  in  the  same  plane  perpendicular  to  the  nozzle 
axis;  however,  they  are  located  at  different  radii  in  this  plane.   The 
tubes  begin  and  end  with  an  inner  diameter  of  0.3  inches.   At  the  noz- 
zle throat  they  form  an  ellipse  with  major  and  minor  axes  of  0.31  in. 
and  0.19  in.,  respectively.   This  elliptical  cross  section  gradually 
increases  until  at  the  bottom  of  the  tubes  the  major  and  minor  axes 
have  values  of  0.76  in.  and  0.70  in.   A  scale  model  of  one  of  these 
tubes  is  shown  in  Fig.  2.   From  the  discharge  of  the  tubes  the  oxidizer 
enters  a  second  annular  chamber  concentric  with  that  into  which  the 
pump  discharges.   In  this  annular  passage  the  flow  is  turned  to  enter 
a  series  of  60  holes  drilled  radially  in  the  inner  shell  of  the  second 
annulus.   The  oxidizer  is  then  discharged  into  the  primary  combustion 
chamber  through  180  injectors.   A  flow  passage,  not  shown  in  Fig.  1, 
directs  the  combustion  gases  to  a  single-stage  axial  turbine  which 
drives  the  pump  previously  mentioned.   Details  of  this  flow  passage  are 
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given  by  Vavra.   Tests  of  two  turbine  designs  have  been  performed,  the 

2  3 

results  of  which  are  reported  by  Vavra  and  Messegee. 

The  portion  of  the  unit  from  the  inlet  of  the  outer  annulus  to  the 
bottom  of  the  rocket  nozzle  is  called  the  oxidizer  discharge.   The  por- 
tion including  the  remainder  of  the  nozzle  tubes  to  the  injectors  is 
called  the  oxidizer  return.   Secondary  combustion  takes  place  after  the 
turbine  discharge  and  ahead  of  the  nozzle  throat. 

A  scale  model  of  these  flow  passages  has  several  modifications  as 
far  as  shape  and  layout  of  the  actual  unit  is  concerned,  but  mostly  re- 
tains flow  passage  dimensions.   The  model,  shown  assembled  in  Fig.  3, 
is  constructed  of  plastic  and  aluminum  and  weighs  136  pounds.   A  sche- 
matic of  the  flow  passages  is  shown  in  Fig.  4.   The  discharge  and  re- 
turn annuli  are  separated  by  7  parts  of  varying  flow  areas  and  shapes. 
Flow  enters  the  model  through  simulated  pump  vanes  which  impart  a  whirl 
component.   The  whirl  is  removed  by  turning  vanes  and  the  flow  is  dis- 
tributed to  part  1,  shown  separately  in  Fig.  5.   Part  1  has  32  equally 
spaced  groups  of  3  holes  with  a  total  flow  area  of  4.712  sq.  in.   Parts 
2  and  3  are  shown  respectively  in  Figs.  6  and  7.   Each  of  these  parts 
has  52  equally  spaced  single  holes  with  a  total  flow  area  of  2.552  sq. 
in.   Parts  4,  5,  and  6  together  with  2  orifice  plates  simulate  the  noz- 
zle tubes.   An  exploded  view  of  this  assembly  is  shown  in  Fig.  8.   The 
first  orifice  plate,  having  holes  of  diameter  d  ■  0.217  in.  and  dia- 
meter ratio  d/D  ■  6  =  0.723  simulates  the  dpwriiflow  ofufcha  noiale 
tubes.   The  second  plate  with  d  ■  0.216  in.  and  Q    ■  0.720  simulates 
the  up- flow  of  the  nozzle  tubes.   Part  7,  s'hown  in  Fig.  9,  has  32 
equally  spaced  triple  holes  with  a  total  flow  area  of  5.574  sq.  in. 
From  part  7  the  flow  passes  into  the  second  annulus  shown  in  Fig.  10, 


18 


where  it  is  turned  to  enter  the  injectors.   Each  of  the  60  injector 
tubes  has  3  holes  labeled  A,  B,  C,  respectively,  from  the  outer  radius 
inward.   The  tubes  are  labeled  clockwise  from  1  to  60  with  the  first 
tube  just  past  the  nine  o'clock  position  when  the  injectors  are  viewed 
from  the  combustion  chamber  (see  Fig.  36). 

The  model  can  be  divided  into  sub-assemblies  for  testing.   Tests 

on  the  oxidizer  discharge  sub-assembly  including  the  first  annulus 

4 
through  part  5,  excluding  the  second  orifice  plate,  were  made  by  Vavra. 

Preliminary  tests  of  the  return  sub-assembly,  part  6  through  the  in- 
jectors including  the  second  orifice  plate,  were  also  made  by  Vavra. 
This  paper  presents  more  complete  tests  of  the  return  sub-assembly,  in 
addition  to  tests  of  the  fully  assembled  unit. 
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SECTION  3 
TEST  INSTALLATION  AND  INSTRUMENTATION 

Available  as  a  source  of  air  is  a  200  H.  P.  centrifugal  compres- 
sor which  delivers  about  4  pounds  mass  of  air  per  second  at  a  pressure 
ratio  of  about  2.   Part  of  this  flow  passes  through  a  cooler;  the  rest 
discharges  through  a  butterfly  valve  to  atmosphere.   Figure  13  is  a 
schematic  of  the  flow  measuring  installation.   The  cooled  air  enters 
plenum  1  then  passes  through  a  4  in.  pipe  discharging  into  plenum  2 
shown  in  Fig.  11.   Flow  measurement  is  accomplished  by  either  a  sharp 
edge  orifice  or  a  flow  nozzle  installed  in  the  4  in.  pipe  (Fig.  12). 
The  design  of  the  orifice  and  the  flow  nozzle  is  in  accordance  with  VDI 
standards  ,  and  the  installation  satisfies  or  exceeds  the  specified  re- 
quirements of  pipe  lengths  upstream  and  downstream  of  it.   In  addition 
a  flow  straightener  is  installed  at  the  discharge  of  plenum  1.   The 
orifice  has  a  diameter  d  ■  3.024  in.  and  a  diameter  ratio  p    ■  0.7429. 
The  nozzle  has  d  ■  3.022  in.  and   ft  ■  0.734.   Initial  tests  conducted 
by  Vavra  used  the  sharp  edge  orifice  to  measure  flow  rate.   Prior  to  all 
tests  performed  by  the  writer  the  flow  nozzle  was  installed. 

Appendix  A  contains  the  derivation  of  the  flow  measuring  equation 

for  the  nozzle.   Similar  information  on  the  orifice  is  given  by  Kelly. 

Plastic  hoses  connect  pressure  taps  with  the  measuring  instruments. 

The  pressure  difference  P  -  P,  is  measured  by  a  Meriam  Vernier  Mano- 

u    d 

meter  filled  with  water.   The  gage  pressure  ahead  of  the  nozzle  is 

measured  by  a  similar  manometer  filled  with  mercury.   The  measuring 

accuracy  of  these  manometers  is  about  +  0.002  in.   The  temperature  t 

—  u 

is  measured  with  an  accuracy  of  +  1  F.  by  an  iron-constantan  thermo- 
couple located  48.5  in.  upstream  of  the  nozzle. 
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Within  plenum  2  there  is  arranged  a  bend  and  a  10  deg.  conical  dif- 
fuser  with  a  series  of  screens  at  the  discharge  to  minimize  flow  insta- 
bilities.  A  second  thermocouple  is  situated  within  plenum  2.   Plenum 
pressure  is  read  on  either  a  mercury  filled  manometer  or  a  water  filled 
manometer  bank.   Three  such  manometer  banks  were  used.   Bank  1  is  a 
250  cm.  40-tube  bank  which  can  be  read  to  an  accuracy  of  +  0.1  cm. 
Bank  2  is  a  175  cm.  60-tube  bank  with  a  reading  accuracy  of  +  0.2  cm. 
Bank  3  is  a  92  in.  20- tube  board  that  can  be  read  with  an  accuracy  of 
+0.1  in. 

The  model  is  attached  to  plenum  2  by  means  of  a  steel  plate  visi- 
ble in  Fig.  11.   Mounted  on  this  plate  is  a  wooden  nozzle  of  circular 
cross  section  to  produce  uniform  entrance  conditions.   Situated  at  the 

throat  of  the  nozzle  is  a  Kiel  probe  measuring  total  pressure  P 

in 
For  tests  with  the  entire  model  assembly  the  conduit  shown  in  Fig.  14 

is  used  to  introduce  air  from  opposite  directions  into  the  first  annu- 

lus.   Static  pressure  taps  measuring  P  are  located  1.5  in.  ahead  of 

the  model  inlet .   In  the  same  plane  are  located  probes  to  measure  P   . 

o 
One  such  probe  is  pictured  in  Fig.  15  and  shown  schematically  in 

Fig.  16.   The  design  allows  rotation  of  the  probe  into  the  flow,  en- 
suring an  accurate  total  pressure  reading. 

Figure  4  shows  the  location  of  the  measuring  planes  within  the 
model.   Pressure  taps  measure  static  or  total  pressure  depending  on 
their  location.   In  each  of  the  planes  Z,  B,  C  and  E  are  arranged 
8  equally  spaced  static  taps.   Planes  Y  and  G  each  have  8  equally 
spaced  total  taps.   Planes  A  and  F  have  4  static  and  4  total  taps, 
spaced  alternately.   Plane  H  has  12  equally  spaced  static  taps.   Plane 
X  contains  no  taps,  but  the  pressure  can  be  accurately  predicted  at 
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this  location.   Total  and  static  pressures  are  measured  at  selected 
injector  nozzles  in  plane  I  by  the  probes  shown  in  Fig.  17.   These 
probes,  having  a  15-deg.  chamfer  at  the  inlet,  were  designed  to  mea- 
sure total  pressure  of  flows  with  yaw  angles  up  to  20  deg.   From 
Vavra's  results  it  was  suspected  that  the  flow  did  not  leave  each  in- 
jector axially.   Therefore  the  new  probe  assembly  was  designed  to  per- 
mit quick  adjustment  to  reach  the  condition  depicted  in  Fig.  18.   Total 
and  static  pressures  at  seven  locations  within  the  second  annulus  were 
obtained  by  inserting  a  removable  probe  which  could  be  oriented  into 
the  direction  of  flow.   All  pressure  taps  were  connected  to  the  mano- 
meter bank  by  plastic  tubing  and  measured  against  atmospheric  or  ple- 
num pressure  depending  on  the  test  configuration. 

Testing  of  the  return  sub-assembly  required  a  modification  of  the 
installation.   An  adapter  fitted  with  screens,  shown  in  Fig.  20,  was 
mounted  on  plenum  2.   Attached  between  the  adapter  and  the  sub- 
assembly was  a  ring  of  52  inlet  nozzles  to  introduce  the  flow  uniformly 
to  the  orifice  ring.   Plane  D,  only  used  with  this  configuration,  has  4 
pressure  taps.   Because  flow  is  essentially  stagnant  within  the  adapter, 

it  is  assumed  that  P„  ■  P   .   A  view  of  this  test  installation,  without 

D    CD 
the  pressure  taps  connected,  is  shown  in  Fig.  21. 

The  inlet  pressure  for  both  configurations  was  remarkably  steady 
and,  except  where  noted,  the  pressures  within  the  model  did  not  fluc- 
tuate more  than  +  0.5  cm.  of  water  even  at  the  highest  flow  rates. 
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SECTION  4 
THEORETICAL  CONSIDERATIONS 

If  a  steady  incompressible,  adiabatic,  viscous  flow  exists  in  a 
pipe,  Bernouillis  equation  between  two  points  on  a  streamline  can  be 

written  as 

2  2 

Vl  V2 

(144)Pl  +  f  J?"  +  e  zl  '  <144>P2  +  ?  2T  +  e  22  +  ^  1-2        <X> 

where : 

2 
P  =  static  pressure   (lb/in  ) 

£  ■  density   (lb  /ft3) 

V  =  velocity  (ft/sec) 

g  -  32.17   (ft/sec2) 

z  ■  height  above  reference   (ft) 

II  »  head  loss   (lb/ft2) 

Since  the  first  three  terms  on  each  side  of  equation  (1)  represent 

the  total  pressure,  the  head  loss  is  equal  to  the  total  pressure  loss 

£P     in  the  pipe. 
1-2 

Experience  shows  that  this  loss  can  be  expressed  by  the  function 

£P     =  f(  a    ,  u   ,  V,  X,   d,  e)  (2) 

1-2      V 

2 
where:   yL(  ■  viscosity  (lbsec/ft  ) 

Jl    -   characteristic  length  (ft) 
d  ■  characteristic  diameter  (ft) 

e  =*  roughness  parameter  based  on  average  variation  of  inside 
radius   (ft) 

Buckingham's  7T  theorem  states  that  the  number  of  independent 
dimensionless  groups  used  to  describe  a  phenomenon  with  n  vari- 
ables is  equal  to  n-r  where  r  is  the  number  of  basic  dimensions.   For 
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the  basic  dimensions  of  force  F,  length  L,  and  time  T,  a  dimensional 
analysis  of  equation  (2)  gives 


AP 


V2     .    jevd    i 


2,  \   %U    '     d  '    d 


e  vz/g        i v 


-  f     1  TT-»  "T>  7  1  (3) 


Equation  (3)  has  four  independent  parameters 

1-2 

=  Euler  number 


v2/g 

^  Vd 


*A 


Reynolds  number  ■  Re 


—  ■  relative  roughness 


—  =»  geometrical  factor 
a 


Usually  equation  (3)  is  expressed  by 


*V»-  ^3  tf  <*> 


t. 


where  <T  is  the  friction  loss  coefficient  that  depends  on  the  relative 

roughness  and  Re.   In  pipe  flow  with  Re  higher  than  about  10   the 

coefficient  XT  remains  essentially  constant.  In  addition,  the  so- 
called  momentum  loss  occurs  if  pipes  have  either  sudden  contractions 

or  expansions.   It  is  primarily  due  to  eddying  turbulence  and  has  been 
found  experimentally  to  satisfy  the  relation 

2 


AP 


h-2-   ^f2T  (5) 


where  /.  is  the  momentum  loss  coefficient  which  depends  primarily  on 
flow  geometry  and  remains  practically  independent  of  Re. 
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In  a  flow  channel  where  frictional  and  momentum  losses  occur,  the 
two  loss  coefficients  may  be  combined  to  give 

where  ~VjT  is  a.  total  loss  coefficient  which  depends  on  Re.   If  Re  is 

fficiently  high,  the  coefficient  y  is  a  function  of  geometry 


su 


alone. 

The  analysis  can  be  extended  to  systems  of  flow  passages  by  con- 
sidering an  equivalent  pipe  for  which  equation  (6)  is  valid.   If  such 

t 
a  system  is  modeled  to  scale,  the  factors  "U/"  and  -?  are  contained  in 

one  loss  coefficient  P  .   For  air  as  a  testing  fluid,  care  must  be 

taken  to  maintain  low  Mach  numbers  M  in  the  model  as  equation  (6)  holds 

for  incompressible  flows  only. 

For  air  at  low  M  and  with  the  velocity  V,  at  point  1  as  a  reference 


Vl  -  -k 

where:   w  =  flow  rate  (lb  /sec) 

m 

A..  ■  area   (ft  ) 
Equation  (6)  may  then  be  written  as 

.  2 


APt 


1-2=*  I  1-2— 2  <8> 

t   Al 

Using  the  equation  of  state  for  a  perfect  gas 

.  ?tl  (144)  (J) 

v     RT  v  ' 

Cl 


25 


where  R  »  53.35 


ft-lb 

lb  °R 

■ 


T   ■  total  temperature  (  R) 
1 


equation  (8)  becomes 


AP. 


1-2 


2gA/   3  1- 


RT. 


.2 
w 


(10) 


The  so-called  equivalent  flow  rate  referred  to  total  conditions  at 
1  is  defined  by 


w 


i   f\ 


IT 


7        (in) 


(ii) 


Equation  (10)  then  becomes 


AP. 


1-2 


-  k 


1-2 


Kl 


(12) 


where  k 


1-2 


fi. 


2/2A, 


(in'4) 


w    will  be  defined  as  the  equivalent  flow  rate  referred  to  static 

* 
conditions.   Unless  otherwise  stated  w  refers  to  total  conditions. 

3, 


Fo 


r  a  liquid  of  specific  weight  o        (lb/ft  )  equation  (6)  becomes 

Li 


AP  -  \ 

1-2 

Regrouping   terms 


?Hf 


.2 


AP. 


1-2 


1-2 


gA2VL 


.2 

w 


(13) 


(14) 
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Determination  of  k1  9  by  air  tests  of  a  model  makes  it  possible  to 
predict  pressure  losses  in  a  unit  for  particular  fluids  at  a  given  flow 
rate* 

Prior  to  the  application  of  this  theory  to  the  ARES  model  it  was 
decided  to  construct  an  arbitrarily  shaped  flow  passage  which  could  be 
tested  with  air  and  water  varying  Re  and  M.   This  passage,  which  is 
shown  in  Fig.  22,  has  both  friction  and  momentum  losses.   It  was  de- 
signed to  operate  in  the  same  Re  range  for  both  fluids.   Its  inlet  and 
discharge  areas  are  equal  to  4  sq.  in.   At  the  inlet  is  fitted  a  con- 
toured nozzle  to  minimize  entrance  losses  into  the  first  section  of  the 
pipe.   Internally  the  model  has  a  sudden  expansion,  a  flow  straightener, 
and  a  sudden  contraction.   The  minimum  flow  area  of  2.23  sq.  in.  occurs 
in  the  flow  straightener. 

Air  was  obtained  from  the  installation  previously  described.   The 
model  was  attached  to  plenum  2,  as  shown  in  Fig.  23.   Inlet  total  pres- 
sure is  measured  by  the  Kiel  probe,  and  inlet  temperature  was  taken  as 

that  measured  in  plenum  2.   Discharge  pressure  P     was  measured  by 

dis 
one  of  the  probes  shown  in  Fig.  17.   Since  conditions  at  the  discharge 

were  found  to  be  nearly  uniform,  except  in  close  proximity  to  the  walls 

of  the  pipe,  the  average  of  five  traverse  points  was  taken  as  P 


dis 


Static  pressure  P..   was  taken  to  be  atmospheric, 
dis 


Values  of  k.   , .  ,  inlet  to  discharge  pressure  loss  coefficient, 
m-dis 

were  established  for  Re  (based  on  d.  )  between  1.58  x  10  and  3.67  x 

10  and  are  plotted  in  Fig.  24  with  the  curve  giving  the  best  fit 

through  the  data  points.   Variation  of  k.  _,.   is  from  -(J. 042  to 

+0.057  k   .   Higher  values  of  Re  could  not  be  reached  due  to  a  lim- 
avg 

iting  manometer  height.   The  highest  Mach  number  M,  obtained  at  the 
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discharge  as  determined  by  the  static  to  total  pressure  ratio,  was 

0.28.   A  slightly  higher  M  was  undoubtably  present  in  the  flow 

straightener. 

Figure  25  shows  the  test  installation  for  the  water  tests  of  this 

model.   At  the  4  in.  discharge  of  a  fire  hydrant  were  available  about 

400  gal.  H-O/min.  at  a  pressure  of  50  psia.   After  steady  flow  was 

reached,  the  discharge  was  deflected  to  fill  a  50  gallon  drum.   By 

weighing  the  drum  before  and  after  each  run  and  knowing  the  elapsed 

time,  the  flow  rate  was  found.   Inlet  static  pressure  P.   was  measured 

in 

just  after  the  entrance  nozzle  by  a  50  psig.  pressure  gage. 

Values  of  k.   ,.   were  calculated  for  Re  between  1.33  x  10  and 
m-dis 

4.88  x  10  for  4  separate  runs.   The  method  of  data  reduction  is  given 

in  Appendix  B,  and  the  results  are  shown  in  Fig.  24.   Except  for  a  few 

points  there  is  excellent  agreement  between  water  data  and  the  curve 

established  by  the  air  test.   A  maximum  deviation  of  10  per  cent  is 

noted  while  most  points  fall  within  2  per  cent  of  the  curve  obtained 

from  the  air  tests.   At  both  extremes  of  the  Re  spectrum  the  pressure 

P.   tended  to  fluctuate,  which  accounts  for  some  of  the  scatter  in 
in 

these  areas.   What  occurs  at  higher  Re  is  hard  to  predict  although  it 
appears  that  k    ,   is  tending  to  approach  some  constant  value. 
Because  of  limitations  imposed  by  the  method  used  to  collect  and  weigh 
water,  the  higher  Re  range  could  not  be  tested.   If  a  flow  measuring 
orifice  were  installed  upstream  of  the  model  inlet  so  that  high  flow 
rates  could  be  measured  accurately,  data  at  high  Re  could  be  obtained. 
These  tests  prove  that  the  use  of  equations  (12)  and  (15)  for  the 
evaluation  of  the  pressure  loss  coefficient  for  liquid  flows  from  air 
tests  is  justified. 
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SECTION  5 
DESCRIPTION  OF  TESTS 

In  addition  to  the  two  configurations  of  the  ARES  model  and  the 
flow  channel  mentioned  previously,  tests  were  made  on  a  full  scale 
model  of  a  nozzle  tube  to  determine  a  representative  loss  coefficient. 
Moreover,  separate  tests  were  made  with  a  geometrically  similar  model 
of  a  tube  section  in  the  ARES  model  from  plane  W  to  plane  X,  simu- 
lating the  nozzle  tubes.   These  tests  are  described  in  the  sequence 
they  were  made.   Table  I  gives  an  index  of  the  runs  conducted.   All  raw 
data  from  these  tests  are  on  file  with  Professor  Vavra  at  the  Pro- 
pulsion Laboratory. 

5. 1  Oxidizer  Return  Sub-assembly 

In  testing  the  oxidizer  return  sub-assembly,  a  total  of  six 
runs  was  carried  out  at  Re  between  0.868  x  10  and  1.343  x  10  referred 
to  conditions  just  upstream  of  the  flow  measuring  orifice.   The  maximum 
Re  attainable  was  restricted  by  the  height  (250  cm)  of  manometer  bank  1. 
As  a  result  of  Vavra' s  preliminary  tests  it  was  determined  that  a  new 
discharge  total  probe  assembly  should  be  designed.   Runs  Rl  and  R2 
were  conducted  without  taking  any  measurements  in  discharge  plane  I. 
The  new  probe  assembly  shown  in  Fig.  17  was  built,  calibrated,  and 
used  first  during  run  R3  where  alignment  and  ease  of  adjustment  were 
investigated.   More  complete  data  in  plane  I  were  taken  during  runs 
R4,  R5  and  R6.   At  each  selected  peripheral  location  the  probes  were 
aligned  individually  with  the  injectors.   Although  the  readings  were 
steady  and  accurate,  the  procedure  was  time  consuming.   It  was  therefore 
decided  to  use  data  from  every  third  injector  tube  only. 
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Pressure  taps  were  made  in  the  annulus  by  drilling  a  1/16  in. 
hole  partially  tapped  with  a  10-32  thread.   If  the  portable  probe  shown 
in  Fig.  19  was  not  inserted,  these  holes  were  plugged  by  screws. 
Because  of  wear  on  the  threads  in  the  plastic,  conditions  in  the  annu- 
lus were  determined  only  for  runs  R3  through  R6. 

Pressures  from  plenum  2  downstream  through  plane  H  were  measured 
on  manometer  bank  1  referenced  to  atmosphere.   Pressure  in  plane  I  was 
measured  on  bank  3  also  referenced  to  atmosphere.   Because  of  the 
large  number  of  tubes  and  the  slight  pressure  fluctuations  noticed  dur- 
ing testing,  Polaroid  pictures  were  taken  of  bank  1  for  more  accurate 
data  recording. 

At  locations  5,6,  7  and  8  in  plane  H  the  pressures  were  observed 
to  fluctuate  considerably,  especially  at  higher  flow  rates.   The  area 
of  this  unsteadiness  is  visible  in  Fig.  10.   For  runs  Rl  and  R2  the 
fluctuation  was  dampened  by  clamping  the  pressure  lines  leading  from 
these  taps  to  the  manometer  bank.   During  runs  R3  through  R6  the  clamps 
were  removed  and  an  average  was  taken  between  the  maximum  and  minimum 
fluctuation.   This  method  produced  the  same  readings  as  with  the  clamps, 
but  the  fluctuations  could  now  be  observed. 

The  general  procedure  for  conducting  each  run  was  as  follows: 

1.  The  micromanometers  and  thermocouples  were  zeroed  and 
atmospheric  pressure  recorded. 

2.  The  desired  pressure  in  plenum  2  was  set  and  the  flow  was 
allowed  to  stabilize. 

3.  Flow  measuring  data  were  recorded  and  a  photograph  of  bank  1 
was  taken. 

4.  Plane  I  and  annulus  data  were  recorded. 

5.  Part  3  was  repeated. 
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Because  of  the  interesting  velocity  distribution  in  the  discharge 
plane  and  the  large  pressure  fluctuations  at  locations  in  plane  H, 
attempts  were  made  to  observe  the  flow  in  the  annulus  visually.  Dif- 
ferent methods  to  introduce  smoke  had  negative  results.   The  first 
attempt  used  a  50  gallon  drum  within  which  a  standard  Navy  MK-5  day 
signal  flare  was  ignited.   This  particular  flare  produces  about  15 
minutes  of  a  dense,  slightly  toxic,  white  smoke.   The  drum  was  con- 
nected directly  to  plenum  2  and  pressurized  so  that  the  smoke  would  go 
directly  into  the  plenum.   At  the  discharge  of  the  model  there  was 
arranged  a  device  to  pipe  the  flow  out  of  the  building.   It  was  found 
that  too  much  smoke  was  produced,  and  the  mixing  in  plenum  2  was  so  com- 
plete that  there  was  no  streamline  definition  in  the  annulus.   The 
flare  produced  enough  heat  to  melt  the  rubber  pressure  seals  in  the 
drum,  allowing  toxic  smoke  to  enter  the  room.   For  these  reasons  this 
attempt  was  abandoned  in  favor  of  introducing  smoke  directly  into  the 
annulus.   Cigars  soaked  in  mineral  oil  and  paraffin  were  ignited  in  a 
small  container  which  was  pressurized.,   Lines  from  this  container  led 
to  eight  1/16  in.  diameter  holes  spaced  around  the  entrance  to  the 
annulus.   The  smoke  partially  condensed  in  the  lines  and  that  part  which 
entered  the  annulus  diffused  quite  rapidly.   Although  streamlines  were 
faintly  discernible,  they  were  of  insufficient  definition  to  observe  the 
flow. 

This  latter  method  left  a  large  residue  of  tar  and  nicotine  in  the 
model.  With  the  air  passing  through  water  was  introduced  into  the 
annulus  to  clean  out  this  residue.   It  became  evident  then  that  water 
would  work  as  a  visualization  fluid.   The  set-up  shown  in  Fig.  21  con- 
sists of  a  manifold  connected  to  a  plastic  liquid  soap  bottle.   Water 
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colored  by  food  dye  was  forced  into  the  annulus  through  the  top  four 
holes  by  squeezing  the  soap  bottle.   This  method  was  very  successful 
and  produced  a  flow  picture  of  photographic  quality.  Motion  pictures 
of  runs  R3  through  R6  were  taken  using  a  Kodak  Cine-Special  II,  16  mm 
camera.   Both  black  and  white  and  color  film  were  used  at  64  frames  per 
second.   Runs  R3,  R4  and  R5  were  photographed  in  black  and  white  using 
TRI-EX,  ASA  120  film  and  a  red  filter.   KODACHROME  II,  ASA  25,  was 
used  on  runs  R5  and  R6. 

5.2  Nozzle  Tubes 

The  nozzle  tubes  were  tested  in  the  configuration  shown  in 
Fig.  26.   Water  was  supplied  by  a  1.0  in.  hose  from  a  faucet  where  the 

pressure  was  higher  than  50  psig.   A  total  of  22  data  points  during  two 

4  4 

runs  was  taken  between  Re  of  1.844  x  10  and  6.874  x  10  ,  referred  to 

inlet  diameter. 

Static  pressure  was  measured  at  the  tube  inlet  by  a  50  psig  pres- 
sure gage  accurate  to  within  +  0„1  psi.   The  flow  discharged  to  atmo- 
sphere-  For  each  data  point  the  inlet  pressure  was  set  and  the  flow 
allowed  to  stabilize.   Water  was  then  collected  in  a  five  gallon  bucket 
during  a  measured  time.   Data  recorded  were  inlet  static  pressure, 
weight  of  water,  and  elapsed  time.   The  maximum  Re  was  obtained  with 
the  water  supply  fully  open. 

5. 3  Geometrically  Similar  Model  of  Plane  W  to  Plane  X 

A  single  tube  from  plane  W  to  plane  X  was  scaled  up  by  a  fac- 
tor of  3.5  and  tested  with  water  by  using  the  set-up  shown  in  Figs.  27 
and  28.   Water  was  supplied  through  a  1  in.  hose  to  a  24  in.  length  of 
6  in.  diameter  pipe,  at  the  bottom  of  which  is  installed  a  contoured 
nozzle  to  ensure  smooth  entry  into  the  1.05  in.  diameter  pipe  section. 
For  runs  Wl  and  W2 ,  the  orifices  A  and  B  had  diameters  of  0.761  and 
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0.757  in.,  respectively.  Prior  to  run  W3  new  orifices,  each  having  a 
diameter  of  0.714  in.,  were  installed.   A  constant  water  level  in  the 
reservoir  was  maintained  at  all  times,  and  the  flow  rate  was  con- 
trolled by  regulating  the  water  supply  and  the  discharge  valve.   The 

Reynolds  number  Re,  referred  to  inlet  diameter,  could  be  changed  from 

4  4 

0.415  x  10  to  4.078  x  10  during  the  tests.   The  runs  were  conducted 

in  the  same  manner  as  those  for  the  nozzle  tubes.   The  pressures  at 
the  6  locations  indicated  in  Fig.  27  were  read  from  a  manometer  board 
referenced  to  atmosphere. 

5.4  Discharge  and  Return  Assembly 

Based  on  the  information  from  these  water  tests  two  new  ori- 
fice plates,  with  holes  of  0.204  in.  diameter,  were  installed  in  the 
ARES  model.   A  total  of  five  runs,  R7  through  Rll,  was  carried  out  at 
Re,  referred  to  the  4  in.  supply  pipe,  between  0.917  x  10  and  1.32  x 
10  .   Pressures  from  plenum  2  downstream  through  plane  B  were  measured 
on  bank  1  referenced  to  plenum  2.   Bank  2,  referenced  to  atmosphere, 
measured  pressures  from  plane  C  through  plane  Ho   Bank  3  was  used  again 
for  plane  I.   The  test  procedures  were  similar  to  those  used  for  the 
return  sub-assembly.   Two  cameras  were  used  to  photograph  manometer 
banks  1  and  2  simultaneously.   No  data  was  taken  within  the  second 
annulus.   Total  pressure  in  plane  Z  was  measured  by  the  removable  probe, 
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SECTION  6 
RESULTS  AND  DISCUSSION 

The  oxidizer  return  sub-assembly  was  tested  separately,  mainly  to 
investigate  the  velocity  distribution  in  the  discharge  plane  and  the 
stability  of  the  flow  inside  the  return  annulus.   To  a  large  extent 
these  criteria  determine  the  feasibility  of  the  flow  passage  design. 
Assuming  that  fuel  is  pumped  into  the  combustion  chamber  at  plane  I  and 
distributed  uniformly  around  the  periphery,  the  importance  of  a  uniform 
oxidizer  distribution  can  be  seen.   The  mixture  of  fuel  and  oxidizer 
will  be  rich  or  lean  depending  on  whether  the  oxidizer  flow  is  defi- 
cient or  excessive.   If  the  oxidizer  discharge  varies  either  with  time 
or  peripherally,  and  if  it  is  not  possible  to  meter  the  fuel  accord- 
ingly, a  non-uniform  combustion  will  result.   This  condition  is  unde- 
sirable for  two  reasons.   First,  hot  spots  will  result  in  the  combus- 
tion chamber  inducing  thermal  stresses  in  the  combustion  liner.   Second, 
and  more  important,  is  the  effect  on  turbine  inlet  conditions.   The 
combustion  chamber  length  is  short,  and  non-uniform  temperatures  would 
have  insufficient  time  to  adjust  prior  to  reaching  the  turbine  inlet. 
A  turbine  blade  then  would  pass  through  a  non-uniform  temperature  field. 
An  off-design  condition  would  result  and  could  have  a  detrimental  effect 
on  turbine  performance.   A  major  portion  of  the  regenerative  cooling  in 
the  ARES  system  occurs  within  the  return  annulus.   The  convective  heat 
transfer  in  this  annulus  is  totally  dependent  on  the  oxidizer  flow 
pattern. 

Since  the  return  sub-assembly  differs  from  the  full  model  only 
upstream  of  plane  E,  it  seems  reasonable  to  assume  that  flow  downstream 


34 


of  this  point  will  represent  the  conditions  that  will  exist  in  the 

fully  assembled  model.   Actually  it  was  found  that  inlet  conditions  at 

plane  E  were  essentially  the  same  in  both  test  configurations. 

Data  reduction  for  runs  Rl  through  R.6  is  contained  in  Appendix  C, 

part  2.   The  equivalent  flow  rate  w   referred  to  conditions  at  the 

flow  measuring  nozzle,  and  the  values  of  Re  just  upstream  of  this  point 

are  listed  in  Table  II.   Appendix  A  gives  the  derivation  of  the  flow 

measuring  equations.   Actual  flow  rate  w  varies  from  0.283  lb  /sec  in 

m 

run  R3  to  0.449  lb  /sec  in  run  R6.   Also  listed  in  Table  II  are  average 
m 

total  pressures  and  equivalent  flow  rates  at  intermediate  planes  D,  F 
and  G  found  according  to  equation  (C2).   Losses  between  planes  will  be 
referred  to  one  of  these  average  inlet  conditions. 

It  is  important  that  the  flow  is  distributed  evenly  around  the 
periphery  as  it  enters  the  model.   In  this  case  the  procedure  of  aver- 
aging conditions  at  each  plane  can  be  justified.   From  water  tests  con- 

4 
ducted  by  Vavra ,   on  a  geometrically  similar  model  of  a  tube  from  plane 

D  to  plane  X,  an  orifice  discharge  coefficient  K  was  obtained.   Vavra 's 

model  scaled  the  orifice  size  by  a  factor  of  3.5  from  0.216  in.  to 

0.757  in.   This  discharge  coefficient,  analagous  to   e<  f  of  equation 

(A3) ,  takes  into  account  some  of  the  difference  between  measured  and 

actual  flow  rate.   Figure  29  shows  K  plotted  as  a  function  of  Re 

referred  to  tube  inlet  conditions.   Values  of  K  in  Fig.  29  appear  some- 

Q 

what  higher  than  those  listed  by  ASME   for  a  sharp  edge  orifice  with 
ft  =  0.7209.   This  is  quite  reasonable  when  the  different  inlet  con- 
ditions are  considered.   The  ASME  data  are  dependent  on  a  minimum 
length  of  straight  pipe  upstream  of  the  orifice.   For  Vavra ' s  tests  an 
inlet  nozzle  was  present  only  ahead  of  the  orifice.   In  later  tests  of 
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the  same  orifice  with  some  length  of  pipe  upstream,  values  of  K  were 

obtained  which  agreed  more  closely  to  those  listed  by  the  ASME. 

The  flow  at  plane  D  was  found  to  be  stagnated,  and  the  pressure 

P  =  P   was  uniform.   At  each  measuring  point  in  plane  E  an  equiva- 

lent  flow  rate  w    was  found  according  to  equations  (C5)  and  (C6) . 

Assuming  that  w   existed  through  each  of  the  52  orifices  the  peri- 

i 
pheral  flow  distribution  can  be  expressed  by  a  factor  F. ,  where 

(52)  w/ 

Fi  "  T-1-  (15) 

WD 

* 
The  total  equivalent  flow  rate  w   is  given  in  Table  II.   Tables  III 

and  IV  list  values  of  F.  in  addition  to  intermediate  quantities  used 
in  the  calculation.  Figure 31  is  a  plot  of  F.  as  a  function  of  peri- 
pheral location  for  Rl  through  R6.   Looking  downstream  the  angle 
is  measured  counterclockwise  with  Q   ■  0  deg.  corresponding  to  the 
three  o'clock  position.   In  each  instance  the  flow  distribution  is 
almost  identical.   The  maximum  deviation  of  this  average  F.  from 
F.  a  1.0  occurs  in  run  R3  where  the  lowest  Re  was  encountered.   It  was 
observed  that  most  of  the  questionable  data  obtained  during  these  tests 
occurred  at  the  lower  Re  where  the  measuring  accuracy  is  lower.   Gen- 
erally each  curve  has  two  peaks  and  two  valleys  indicating  a  deficiency 
of  flow  at  0  =  110  deg.  and  300  dego  and  an  excess  at  9=0  deg.  and 
225  deg.   This  is  a  relatively  symmetric  condition  with  the  maxima 
appearing  about  180  deg.  apart.   Because  of  this  somewhat  regular 
occurrence  the  values  of  F.  at  Q   ■  90  deg.  for  runs  R3  through  R6  are 
questionable.   Not  be  be  overlooked,  however,  is  the  fact  that  the 
maximum  deviation  from  F.  *  1.0  is  only  4  per  cent.   This  fact  tends  to 
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establish  the  accuracy  of  Fig.  29  and  to  substantiate  that  flow  enters 
the  oxidizer  return  sub-assembly  uniformly. 

Figure  30,  also  plotted  from  Vavra's  water  tests,  shows  j    , 
the  loss  in  total  pressure  from  plane  D  to  plane  X  divided  by  the  static 
pressure  change  across  the  orifice,  as  a  function  of  Re.  The  total  pres- 
sure P   at  plane  X  is  found  according  to  equation  (C8).   With  this 

X 
pressure  and  the  inlet  conditions  at  plane  D,  the  coefficient  k    was 

calculated  for  each  peripheral  location.   An  average  value  k    for 
each  run  was  found  which  is  shown  in  Tables  III  and  IV.   It  should  be 
noted  that  at  any  peripheral  location  the  value  of  k^   differs  from 
k^  Y  by  not  more  than  five  per  cent. 

If,  as  was  hypothesized,  the  losses  in  this  section  are  indepen- 
dent of  Re,  then  k^   should  not  vary  from  point  to  point.   This  may  be 

shown  by  referring  to  equation  (12)  and  recognizing  that  both  AP 

*2  - 

and  w_    depend  essentially  on  AP   .   The  distribution  of  k^  v/k_  v 
D,  D-fci  l)-X   D-X 

follows  the  same  pattern  as  that  established  for  F . .   This  is  shown  in 

* 
Fig.  31  for  Rl  and  R2  and  holds  for  the  other  runs  also.   Since  w   is 

Di 

proportional  to  Re  in  accordance  with  equation  (C6) ,  the  loss  coeffi- 
cient is  also  proportional  to  Re. 

Experimentally  obtained  values  of  the  friction  factor  proportional 
to  k  as  a  function  of  Re  and  relative  roughness  are  shown  in  Fig.  32 

for  two  different  conditions  of  flow  in  pipes.   These  curves  were  taken 

9 
from  SchliShting.    For  the  artificially  roughened  pipe  of  Fig.  32(a) 

the  friction  factor  increases  with  Re  to  approach  a  constant  value, 

whereas  for  the  commercial  pipe  of  Fig.  32(b)  the  factor  decreases  with 

Re  toward  a  constant  value.   In  the  former  case  the  pipes  tested  were 

lined  with  protuberances  which  penetrated  completely  through  the 
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laminar  sub-layer  of  the  boundary  layer.   A  large  part  of  the  friction 
is  then  due  to  a  form  drag  which  can  be  thought  of  as  a  momentum  loss. 
In  the  second  case  skin- friction  alone  accounts  for  the  friction  fac- 
tor.  As  relative  roughness  is  increased  in  both  cases  the  friction 
factor  reaches  a  constant  value  at  a  lower  Re.   Realizing  that  in  the 
present  model  test  each  tube  from  plane  D  to  plane  X  is  machined  to  at 
least  the  equivalent  of  commercial  roughness  and  that  k^   increases 
slightly  with  Re,  it  can  be  concluded  that  the  momentum  loss  created  by 
the  orifice  predominates.   Tables  III  and  IV  show  values  of  Re  corres- 
ponding to  the  general  transition  range  of  Fig.  32.   The  loss  coeffi- 
cient k^   would  become  constant  if  either  more  momentum  loss  were  in- 
troduced or  if  Re  were  increased. 

Incompressible  flow  is  a  basic  assumption  of  the  theory  proposed 
for  this  analysis.   It  seems  reasonable  to  assume  that  the  highest  Mach 
number  M  in  the  model  exists  at  the  minimum  flow  area  which  occurs  at 
the  vena-contracta  downstream  of  the  orifice  plate.  At  plane  E  an 
average  area  restriction  factor  Cp       was  found  using  equations (C9) 

through  (C13).   Also  calculated  were  average  values  of  P  ,  P   and  M  . 

E 
Table  V  lists  these  values  for  runs  Rl  through  R6.   Applying  the  values 

of  0   to  a  flow  area  at  plane  E  of  3.676  sq.  in.  gives  an  effective 

flow  area  varying  from  2.081  sq.  in.  to  2.161  sq.  in.   The  next 

largest  area  in  the  model  is  the  3.138  sq.  in.  area  at  plane  I,  where 

an  average  M  was  calculated  and  is  listed  in  Table  V.   For  all  runs 

the  highest  M  in  the  model,  occurs  at  plane  E,  reaching  a  maximum  value 

of  0.34  in  run  R6.   Hence  it  is  permissible  to  assume  that  the  flow  is 

essentially  incompressible. 
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Table  V  also  lists  values  of  average  loss  coefficients  from  planes 
D  to  F,  F  to  G,  G  to  I  and  X  to  F.   The  latter  loss  is  given  by 


:-F  "  ^-F  "  T>-2 


(16) 


From  Fig.  30  are  obtained  values  of  \   which  can  be  used  to  pre- 
dict P   ,  and  hence  k    ,  whereas  the  values  k.^   are  calculated  from 

4m 

test  data.   The  magnitude  of  k    is  small  because  of  the  small  dis- 

X-F 

tance  between  the  two  planes.   It  can  be  said,  therefore,  that  the 
method  used  to  establish  Fig.  30  is  a  valid  one. 

The  peripheral  distribution  of  k^   /k^   is  shown  in  Table  VI. 
Except  for  a  few  isolated  values,  kT)_„  is  contained  within  2  per  cent 
of  kn_F>   Referring  to  the  previous  discussion  concerning  the  rela- 
tionship between  equivalent  flow  rate  and  loss  coefficient,  it  can  be 
assumed  that  the  flow  is  distributed  uniformly  at  plane  F. 

The  situation  is  different  in  plane  G.   Figure  33  shows  the  per- 
ipheral distribution  of  k   /k   .   Since  the  pressure  tap  at  Q  3  0  deg, 

F-G  F~G 

was  damaged,  a  curve  was  faired  in  to  give  an  estimated  value  at  this 
location.   Except  for  runs  R3,  at  low  Re,  a  definite  pattern  can  be 
observed.   A  minimum  loss  coefficient  occurs  at  ft"  180  deg.   This  in- 
dicates an  irregularity  of  the  flow  at  this  peripheral  location.   This 
situation  can  be  explained  from  the  geometry  of  the  flow  passages. 
Figure  34  shows  the  intersection  of  the  two  cylinders  with  the  posi- 
tions of  Q   in  plane  G  indicated.   At  ©  =  90  deg.  and  270  deg.  the 
flow  has  an  easier  path  to  follow,  therefore  one  would  expect  a  higher 
flow  rate  at  these  points  and  hence  a  higher  value  of  loss  coefficient. 
However,  at  &   =  0  deg.  the  flow  enters  the  return  annulus  on  top  of 
the  holes  in  plane  H.   It  will  be  shown  than  the  primary  injectors  at 
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this  location  pass  a  higher  percentage  of  the  discharge  flow.   It  is 

entirely  reasonable  then  that  the  flow  is  deficient  as  it  enters  the 

annulus  from  &    =   135  deg.  to  Q   =  225  deg. 

Table  VII  lists  values  of  loss  coefficient  k    for  each  location 

G-H 

in  plane  H  and  k    for  runs  Rl  through  R6.   Also  shown  in  Table  VII 
G~H 

are  values  of  L   for  runs  R3  through  R6,  found  by 

Vi '  Vi "  Vh  <17> 

Shown  in  Fig.  35  are  values  of  average  loss  coefficients  between 
different  stations  of  the  model  as  a  function  of  Re.   The  increase  with 
Re  of  L    has  already  been  explained,  however,  the  maximum  and  minimum 
values  are  only  92  and  104  per  cent  of  the  average  value,  respectively. 
The  highest  value  which  k^   approaches  at  increased  Re  is  probably 
around  0.125.   If  loss  coefficients  in  the  low  Reynolds  number  range 
are  not  considered,  the  average  losses  from  plane  X  throughout  the  re- 
mainder of  the  model  are  remarkably  constant. 

Velocity  distribution  in  the  discharge  plane  is  found  according 
to  equations  (C15)  and  (C16).   Figure  36  shows  a  drawing  of  plane  I 
from  the  discharge  side  and  the  relative  position  of  the  pressure  taps 
in  plane  H.   Table  VIII  lists  values  of  V . /V  for  each  of  the  surveyed 
locations  for  runs  R3  through  R6.   An  average  of  positions  A,  B  and  C 
was  calculated  at  each  location  1  through  57.   This  average  velocity 
distribution  is  plotted  in  Figs.  37  through  40.   The  velocity  distribu- 
tion and  hence  the  flow  distribution  is  essentially  independent  of  Re. 
The  excessive  flow  from  locations  54  through  10  was  anticipated  from 
the  conditions  existing  at  0  =  0  deg.  in  plane  G-   The  erratic  con- 
dition from  locations  24  through  36  will  be  explained  later.   It  is  not 
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understood  why  the  flow  is  deficient  around  location  48.   Data  at  the 
seven  locations  within  the  annulus  obtained  from  the  removable  probe 
were  entirely  unreliable.   The  model  had  been  dropped  causing  damage  to 
the  plastic  threads  and  making  it  virtually  impossible  to  prevent 
leakage. 

The  flow  within  the  return  annulus  can  best  be  described  by  re- 
ferring to  the  results  of  the  flow  visualization  efforts  previously  de- 
scribed.  It  has  been  shown  that  the  flow  in  plane  G  is  symmetrical 
about  the  Q    =  0  deg.  and  0  =  180  deg.  axis.   Therefore,  the  annulus 
is  naturally  split  into  a  top  and  bottom  half  in  each  of  which  the  same 
flow  pattern  should  exist.   Water  introduced  in  the  lower  half  had  a 
tendency  to  pool  at  the  lowest  point  because  of  gravity  effects,  whereas 
all  the  water  introduced  in  the  top  half  was  carried  out  of  the  model. 
For  these  reasons  and  because  of  problems  associated  with  photographing 
the  flow  conditions,  it  was  decided  to  investigate  only  the  upper  half 
of  the  annulus.   Figure  41  is  a  drawing  of  that  half  annulus,  developed 
into  a  plane,  showing  the  locations  where  water  was  introduced.   In 
addition,  data  obtained  by  Vavra  are  shown  that  give  flow  directions 
at  locations  within  the  annulus.   These  flow  directions  are  represented 
by  arrows  in  Fig.  41.   Figure  42(a)  through  (h)  and  Fig.  43(a)  through 
(h)  are  photographs  taken  of  runs  R4  and  R5,  respectively.   The  upper 
half  of  the  annulus  is  shown  as  the  upper  3/4  of  each  picture.   These 
pictures  were  made  from  single  frames  of  16mm  movie  film  taken  at  64 
frames  per  second.   Each  sequence  spans  one  second  of  operation  with  a 
time  between  each  pattern  of  0.125  sec.   Figure  42  is  typical  of  the 
flow  pattern  at  Re  ■  1.028  x  10  and  lower> whereas  Fig.  43  is  repre- 
sentative of  the  conditions  at  Re  =  1.177  x  10  and  higher.   Several 


41 


similarities  exist  at  all  Re.   Most  apparent  is  the  extreme  instability 
in  the  region  of  plane  H  where  the  flows  from  the  top  and  bottom  halves 
of  the  annulus  mix.   This  region  is  exactly  that  in  which  the  pressure 
taps  H5,  H6,  H7  and  H8  are  located. 

It  should  be  noted  that  this  is  the  only  region  in  the  model  where 
the  pressures  were  unsteady.   Several  vortices  are  present  in  this  gen- 
eral area.   The  pattern  is  apparently  being  triggered  by  a  huge  separa- 
tion from  the  turning  vane  near  H5.   Although  not  entirely  obvious  from 
the  sequence  of  pictures,  it  was  observed  that  these  vortices  went 
through  a  periodic  variation  in  both  intensity  and  location.   At  best, 
it  can  be  said  that  the  flow  in  this  area  is  thoroughly  confused  and 
that  this  situation  is  most  likely  the  cause  of  the  cyclic  velocity  dis- 
tribution in  plane  I  which  is  evident  from  Figs.  37  through  40.   It  has 
been  noted  also  that  a  deficiency  of  flow  exists  at  the  far  right  of 
the  annulus.   This  portion  of  the  annulus  is  supplied  by  the  Q   ■  180 
deg.  location  of  plane  G  where  the  deficiency  of  flow  has  already  been 
observed.   The  most  striking  difference  between  Fig.  42  and  Fig.  43  is 
the  appearance  in  the  latter  of  a  huge  vortex  in  the  upper  right  hand 
portion  of  the  annulus.   It  is  believed  that  this  vortex  is  caused  by  a 
flow  separation  from  the  turning  vane  immediately  upstream  of  it.   A 
100  ft.  roll  of  color  film  taken  of  runs  R5  and  R6  was  left  intact,  and 
is  on  file  with  Professor  Vavra. 

From  the  above  discussion  and  the  flow  visualization  studies,  it 
is  possible  to  make  a  qualitative  design  evaluation  of  this  portion  of 
the  ARES  assembly.   The  test  data  at  the  higher  Re  will  best  simulate 
actual  operation,  although  this  distinction  appears  necessary  only  when 
considering  conditions  within  the  annulus.   Heat  transfer  will  most 
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certainly  be  hampered  where  vortices  exist,  since  the  same  flow  parti- 
cles will  tend  to  circulate  in  one  spot  instead  of  passing  through  the 
model.   Hot  spots  could  result  in  these  regions.   The  deficiency  of 
flow  in  the  region  farthest  from  plane  H  would  have  the  same  effect  on 
heat  transfer,  only  to  a  lesser  extent.   Elimination  of  these  flow  in- 
stabilities, in  addition  to  improving  heat  transfer  characteristics, 
would  serve  to  reduce  the  losses  from  plane  G  to  H  and  to  distribute 
the  discharge  velocity  more  uniformly.   One  method  of  accomplishing 
this  would  be  to  redesign  the  turning  vanes,  making  their  shape  conform 
more  to  airfoil  cross  sections  cambered  in  the  direction  in  which  the 
flow  is  to  be  turned.   For  stress  considerations,  the  contact  surface 
area  with  the  outer  annulus  would  have  to  be  maintained.   Another 
scheme  would  be  to  replace  the  relatively  large  number  of  turning  vanes 
with  a  few  long  flow  dividers.   It  is  felt  that  a  large  portion  of  the 
instability  occurs  as  a  result  of  the  mixing  of  flows  from  each  half  of 
the  annulus,  and  for  this  reason  the  arranging  of  a  flow  splitter  is 
recommended. 

The  loss  coefficients  calculated  during  the  tests  of  the  return 
sub-assembly  were  duplicated  when  the  full  model  was  tested  and  will  be 
discussed  later  in  this  report.   Before  proceeding  with  the  tests  of 
the  entire  unit,  it  was  decided  to  determine  how  closely  the  pressure 
drop  in  the  model,  from  plane  W  to  plane  X,  simulated  that  in  the 
actual  nozzle  tubes.   Results  of  tests  on  the  nozzle  tube  are  listed  in 
Table  IX.   The  method  of  data  reduction  is  shown  in  Appendix  C,  part  3. 

Loss  coefficient  P  .   ,.   is  plotted  as  a  function  of  Re  in  Fig.  44. 
}  m-dis 

Fluctuation  of  the  pressure  gage  was  as  much  as  +  0.5  psig.   At  low 
values  of  Re,  where  pressures  were  of  the  order  of  3.5  psig.,  this 
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fluctuation  could  cause  the  computed  "P   to  vary  between  11.1  and  8.5. 
Therefore,  the  dashed  portion  of  the  curve  of  Fig.  44  is  questionable. 
The  decrease  of  P  with  Re  is  quite  evident,  and  by  applying  the 
rationale   previously  described,  it  seems  that  skin  friction  accounts 
for  the  major  portion  of  the  loss.   Because  of  the  scale  chosen  for  F 
in  Fig.  44  its  variation  appears  magnified.   Actually  £  only  deviates 
+  6  per  cent  from  an  average  value  of  10.3.   The  range  of  Re  encountered 
falls  well  within  the  transition  range  for  pipes  as  shown  in  Fig.  32. 

Data  from  tests  on  the  geometrically  similar  section  of  plane  W  to 
plane  X  were  reduced  by  the  method  outlined  in  Appendix  C,  part  4. 
These  tests  were  conducted  similarly  to  those  run  by  Vavra  for  which  the 
validity  has  been  established.   Table  X  lists  values  of  t  -  ,,  h  -  , 
and  the  orifice  discharge  coefficients,  in  addition  to  other  quantities 
of  interest.   Points  1  and  6  refer  to  planes  W  and  X  respectively. 
Figure  45  shows  J  , _,  and  F  ,  ,  for  runs  Wl  and  W2 ,  plotted  as  a  func- 
tion of  Re.   The  loss  coefficient  T  -    ,  differs  by  almost  a  factor  of 

2  from  r  .   ,.   of  the  nozzle  tubes.   This  condition  is  unsatisfactory 
)  m-dis  J 

because  the  largest  loss  in  the  entire  assembly  occurs  in  these  tubes 
and  the  water  model  with  orifice  plates  of  d  =•  0.757  in.  and  0.761  in. 
gave  half  of  this  loss  only.   It  is  also  interesting  to  compare  the 
values  of  K  (0.757)  with  those  obtained  by  Vavra.   The  values  differed 
on  an  average  by  almost  6  to  7  per  cent.   The  difference  in  the  dis- 
charge coefficient  is  due  to  the  influence  of  the  inlet  conditions. 

Table  XI  also  lists  the  results  of  run  W3  with  orifice  diameters 
of  0.714  in.,  whereas  Fig.  46  shows  F  .  ,  and  f  .  ,  as  functions  of 
Re.   It  appears  as  if  this  orifice  diameter  more  nearly  simulates  con- 
ditions in  the  nozzle  tube.   From  the  test  of  the  fully  assembled  model, 
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values  of  T     were  computed  as  outlined  in  Appendix  C,  part  5, 
equations  (C32)  to  (C36).   Figure  47  shows  these  values  together  with 
the  curve  obtained  from  the  water  tests  of  the  nozzle  tubes.   The  dif- 
ference in  inlet  conditions  to  plane  W  when  the  section  is  installed  in 

the  model  could  account  for  the  differences  between  f  T.  _  and   ■£,,-. 

j  VJ-X      }  1-b 

Apparently  though  the  orifice  size  of  0.714  in.  does  the  best  job  of 
simulation. 

Tests  of  the  complete  assembly  of  Fig.  4  were  conducted  mainly  to 
establish  values  of  loss  coefficients  without  making  pressure  or  veloc- 
ity surveys.   The  methods  of  data  reduction  are  similar  to  those  used 
for  the  return  sub-assembly  and  are  contained  in  Appendix  C,  part  5. 
Table  XII  gives  the  data  for  the  flow  rate  determined,  together  with 
average  total  pressures  and  equivalent  flow  rates  in  planes  Z,  A,  F 
and  G.   Data  at  plane  Y  indicated  a  lower  total  pressure  than  that 
occurring  in  plane  A.   This  impossible  situation  must  be  due  to  inac- 
curate measurement  in  plane  Y.   These  data  were  therefore  disregarded. 
Average  loss  coefficients  from  planes  Z  to  A,  A  to  F,  F  to  G,  G  to  H, 
H  to  I  and  W  to  X  are  listed  in  Table  XII  for  runs  R7  through  Rll. 

The  range  of  Re  encountered  in  the  4  inch  supply  pipe,  and  in  the 
model  at  plane  E,  was  the  same  as  for  the  return  sub-assembly  alone  as 
was  the  value  of  discharge  Mach  number  M  .   Arranging  the  discharge 
sub-assembly  ahead  of  the  return  sub-assembly  does  not  disturb  the 
peripheral  velocity  distribution  that  was  observed  in  the  return  sub- 
assembly alone.   Table  VI  lists  values  of  k    „/k.   _  which  differ 

A  -  F  A  -  F 

from  a  value  of  1.0  by  a  maximum  of  0.005.   The  distribution  of  loss 
coefficients  from  plane  F  to  plane  G  is  shown  in  Fig.  48.   It  is  almost 
identical  to  that  from  runs  Rl  through  R6.   Velocity  distribution  in 
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plane  I  was  also  found  to  agree  precisely  with  that  determined  pre- 
viously.  From  the  foregoing  it  is  reasonable  to  conclude  that  the  dis- 
cussion of  results  for  the  return  sub-assembly  is  applicable  to  the 
entire  model. 

Figure  49  shows  the  average  loss  coefficients  as  a  function  of 
Re.   The  scatter  of  L   occurs  because  the  total  conditions  at  plane 
A  were  measured  at  only  one  location,  hence  a  true  average  could  not  be 
obtained.   However,  this  condition  is  not  serious  because  this  loss  is 
only  a  small  percentage  of  the  total  loss.   The  total  loss  coefficient 
k,    is  obtained  by  adding  the  coefficients  between  intermediate  planes. 
If  only  the  conditions  in  planes  Z  and  I  are  considered,  a  value  k 
is  obtained  which  differs  from  k^   .   As  Re  is  increased  k^ _T/k7   ' 
goes  from  0.974  to  0.948.   Evidently  there  exists  a  slight  effect  of 
compressibility.   The  difference  in  the  two  loss  coefficients  may  be 
shown  mathematically.   Assume,  for  example,  that  the  loss  coefficient 
between  planes  Z  and  F  is  to  be  determined.   The  overall  coefficient 
may  be  expressed  as 


kz- 


-  *2  - 
P. 


w, 


(Pt  -  Pt  )  +(Pt  "  Pt  ) 
Z     CA      A     F 


(18) 


Using  the  relation 


—  *  —     —  *  — 

WZ   Pt  "  WA  Pt 
L        h  A    CA 


equation  (17)  becomes 


(19) 


_   )   _ 


Z  T 


"z-r "  "*-a  *  r1   k 


A-F 


(20) 
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Since  total  pressure  must  decrease  from  plane  to  plane,  the  pressure 

ratio  P  /P   must  be  larger  than  unity.   Hence,  equation  (20)  shows 

-   ,  A 

that  k^-'  wiH  always  be  larger  than  k^   .   The  total  loss  coefficient 

from  plane  Z  to  plane  I  can  be  expressed  by 

Vl   -  *Z-A  +  Vf  +  Vg  +  Vh  +  W  (2D 

where 

P 

P 


-    '      h  - 

kA-F=~  Vf  (etC'> 

fcA 


Equation  (21)  was  used  to  compute  k^   for  runs  R6  through  Rll.   The 
values  were  found  to  be  identical  to  those  obtained  from 


(22) 


Consider  the  section  from  plane  G  to  plane  H  only.   Figures  35 

and  49  show  that  the  loss  coefficient  as  a  function  of  Re,  if  based 

only  on  inlet  conditions  at  plane  G,  is  the  same  for  both  the  return 

sub-assembly  and  the  fully  assembled  model.   If  the  method  of  equation 

(21)  were  applied,  it  can  be  seen  that  k    would  increase  as  more 

G~H 

sections  were  added  to  the  model  upstream,  because  for  the  same  inlet 

conditions  at  plane  G,  namely,  P   and  w   ,  the  pressure  ratio  P   /P 

fcG      G  h      fcG 

would  increase.   This  method  actually  refers  all  loss  coefficients 

back  to  the  initial  inlet  conditions  and  assumes  incompressible  flow 

from,  say,  plane  Z  to  plane  H.   Using  k    to  represent  the  loss. 

G_H 

effectively  cancels  any  errors  introduced  by  the  incompressible  assump- 
tion for  flows  upstream.   Only  in  that  portionLaf  the  model  from 
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G  to  H  is  incompressible  flow  assumed.   For  this  reason,  it  was  decided 

to  use  the  above-mentioned  simple  method  for  computing  total  loss 

coefficients.   As  an  extension  of  this  reasoning,  the  writer  believes 

that  more  accurate  estimates  of  loss  coefficient  could  be  obtained  by 

adding  more  intermediate  measuring  planes.   It  is  questionable  whether 

this  complication  is  justified  since  the  maximum  difference  between 

k_   and  k_   is  only  5.2  per  cent. 

From  Fig.  49  the  coefficient  k^   appears  to  approach  a  value  of 

about  0.62  of  which  66  per  cent  is  from  plane  W  to  plane  X.   As  an 

example  it  will  be  assumed  that  the  flow  rate  of  the  oxidizer  w  is 

100  lb  /sec  and  that  its  specific  weight  ■  _  is  75  lb/cu.  ft.   The 
m  L 

total  pressure  drop  AP     in  the  oxidizer  passage  from  the  pump  dis- 

Z-I 
charge  to  the  entrance  into  the  primary  combustion  chamber  is  then  from 

equation  (14) 


-%'_x  "  V.  T^f   <lb/i"2>  <"> 


—  -4 

where:  V-        ■  0.62   in 

g    -  32.17  ft/sec2 


or 


AP     -  0.62  (100)2  (144)   ,  37Q  lb/in2  (U) 

Vl    (32.17)   (75) 
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SECTION  7 
CONCLUSIONS 

The  average  loss  coefficient  from  a  position  just  downstream  of 
the  oxidizer  pump  to  the  discharge  of  the  primary  injectors  is  0.62. 
This  experimentally  determined  coefficient  should  remain  nearly  con- 
stant near  the  specified  operating  range  since  the  anticipated  Reynolds 
number  is  sufficiently  high.   The  model  section  from  plane  W  to  plane  X 
with  2  orifice  plates,  having  holes  of  d  ■  0.204  in.,  represents  the 
actual  nozzle  tubes.   The  pressure  loss  coefficient  found  in  the  air 
test  of  this  model  section  agreed  closely  with  the  coefficient  deter- 
mined from  the  water  test  of  the  tubes.   The  peripheral  distribution 
of  flow  within  the  model  is  very  nearly  uniform  until  it  enters  the 
annulus.   Here  the  flow  becomes  deficient  in  a  45  deg.  sector  centered 
on  0  ■  180  deg.  in  plane  G.   This  causes  a  deficiency  of  flow  in  that 
portion  of  the  return  annulus  farthest  from  plane  I. 

Because  of  compressibility  effects  a  difference  exists  between 
loss  coefficients  determined  from  inlet  and  discharge  conditions  only, 
and  those  obtained  by  a  summation  of  the  intermediate  loss  coefficients 
between  adjacent  planes.   The  flow  channel  used  to  check  the  basic 
theory  by  water  and  air  tests  was  neither  designed  nor  instrumented  with 
this  difference  in  mind.   It  would  be  of  interest  to  construct  a  more 
complicated  system  of  flow  passages  for  this  purpose  and  to  investigate 
the  effect  of  compressibility  more  fully.   It  is  believed  by  the  writer 
that  more  intermediate  planes  where  exact  conditions  can  be  established 
will  lead  to  more  accurate  estimations  of  loss  coefficient. 

Due  to  the  fact  that  the  major  portion  of  the  pressure  drop  is 
caused  by  momentum  losses,  it  would  be  inconsistent  to  apply  standard 
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theories  that  relate  heat  transfer  to  pressure  drop  with  the  objective 
of  evaluating  the  heat  transfer  coefficients  at  different  locations  in 
the  flow  passages.   However,  the  vortices  within  the  return  annulus,  as 
evidenced  by  the  photographs  of  the  existing  flow  pattern,  will  cer- 
tainly be  detrimental  to  heat  transfer  in  this  region.   These  vortices 
could  be  largely  eliminated  by  changing  the  design  of  the  turning  vanes 
and  installing  a  flow  divider  where  mixing  now  occurs  between  flows 
from  the  top  and  bottom' of  the  annulus. 

Model  tests  of  this  type  provide  a  natural  bridge  between  the 
completion  of  a  design  and  the  building  of  a  prototype.   It  is  estima- 
ted that  a  cost  factor  of  20  to  1  exists  between  the  actual  hardware  of 
the  ARES  unit  and  its  corresponding  flow  model.   Prediction  of  the 
pressure  loss  coefficients  of  complicated  passages  analytically  is,  at 
present,  an  impossible  task.   Except  for  the  slight  effects  of  Reynolds 
number,  the  experimental  method  of  predicting  pressure  loss  coefficients 
presented  in  this  paper  is  both  valid  and  accurate. 
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TABLE  I 
INDEX  OF  RUNS 

RUN  NUMBER  TEST  DESCRIPTION 

RMl  Air  test  of  flow  device 

RMWl,  RMW2,  Water  test  of  flow  device 

RMW3,  RMW4 

Rl,  R2,  R3,  Air  test  of  oxidizer  return  sub-assembly 

R4,  R5,  R6 

Wl,  W2  Water  test  of  geometrically  similar  model 

plane  W  to  plane  X  orifice  plates  of 
d  =  0.761  in.  and  d  =  0.757  in. 

W3  Same  as  Wl  and  W2  with  orifice  plates  both 

of  d  =  0.714  in. 

WVl,  WV2,  WV3  Water  tests  of  geometrically  similar  model 

plane  W  to  plane  X 

R7,  R8,  R9,  RIO,  Rll       Air  tests  of  full  ARES  assembly 
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TABLE  II 
OXIDIZER  RETURN  SUB-ASSEMBLY,  INLET  ORIFICE 
DATA  AND  EQUIVALENT  FLOW  RATES 

RUN 
AVERAGE  OF  ALL 
DATA  POINTS 

Rl 

R2 

R3 

R4 

R5 

R6 

P  „ 
a  ton 

(lb/in2) 

14.768 

14.768 

14.778 

14.778 

14.792 

14.734 

P 

u 

(lb/in2) 

17.758 

17.550 

16.094 

16.746 

17.437 

18.147 

A  h 
o 

(in  H20) 

1.998 

1.887 

0.971 

1.380 

1.776 

2.221 

T 
u 

(°R) 

539.0 

539.7 

540.5 

541.0 

552.0 

551.5 

plen 

(°R) 

534.7 

535.3 

536.0 

542.0 

550.5 

549.5 

Re(10" 

5) 
'd  -  4.07" 

1.307 

1.262 

0.868 

1.028 

1.177 

1.343 

* 
W 
u 

(in2) 

.71785 

.70176 

.52655 

.61562 

.68344 

.74901 

PD 

(lb/in2) 

17.708 

17.490 

16.058 

16.699 

17.367 

18.061 

\ 

(lb/in2) 

16.510 

16.385 

15.546 

15.926 

16.323 

16.704 

\ 

(lb/in2) 

16.227 

16.122 

15.427 

15.738 

16.073 

16.383 

* 

W 

D 

(in2) 

.71732 

.70137 

.52558 

.61799 

.68534 

.75132 

—  * 
WF 

(in2) 

.76914 

. 74870 

.54291 

.64797 

.72918 

.81237 

_  * 

(in2) 

.78254 

. 76088 

. 54708 

.65574 

. 74051 

.82829 
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TABLE  III 
FLOW  DISTRIBUTION  THROUGH  ORIFICE  PLATE, 
LOSS  COEFFICIENT  PLANE  D  TO  PLANE  X 

RUNS  Rl,  R2,  R3 

Location 

1 

2 

3 

4 

5 

6 

7 

8 

Rl      w  *  =  .71732,  K  -  .79,  P.   -  P_  -  .  50  &   P_  _  =  1161.4 
u                                                           C      D                              D-E 

P  -P 
D  E 

(Cm  HgO) 

163.10 

159.40 

155.05 

151.80 

159.50 

163.10 

159.10 

158.45 

Re(l0~4) 
d  =  0.3 

3.503 

3.468 

3.426 

3.391 

3.468 

3.503 

3.466 

3.458 

* 
w      2 
\   (^  ) 

.01417 

.01403 

.01386 

.01371 

.01403 

.01417 

.01402 

.01399 

F. 

l 

1.027 

1.016 

1.005 

0.994 

1.005 

1.027 

1.016 

1.014 

Vx  (in"4) 

.1273 

.1244 

.1210 
avg. 

.1185 
.1239 

.1245 

.1273 

.1242 

.1237 

R2      wD  -  .70137,  K  =  .79,  Pt  =  PD  -  .50  A  P    =  1152.13 

p  -p 

D  E 

(CmH20) 

151.25 

147.80 

143.85 

141.70 

147.65 

151.60 

147.60 

147.33 

Re(IO"4) 
d  =  0.3 

3.363 

3.328 

3.286 

3.264 

3.326 

3.366 

3.326 

3.324 

*     9 
WDt  (ir/) 

.01377 

.01363 

.01346 

.01337 

.01362 

.01378 

.01362 

.01361 

Fi 

1.021 

1.010 

0.998 

0.991 

1.009 

1.022 

1.009 

1.009 

kD-X(ln:4) 

.1250 

.1221 

.1188 
avg. 

.1171 
.1217 

.1220 

.1252 

.1219 

.1217 

R3      wD  =  .52558,  K  -  .795,  Pfc  -  PD  -  .505  A  P    =  1095.29 

P  -P 
D  E 

(CmHjO) 

69.70 

65.50 

68.95 

64.75 

66.85 

71.45 

67.55 

67.35 

Re  (I0"4) 
d=0.3 

2.249 

2.187 

2.238 

2.173 

2.207 

2.276 

2.218 

2.214 

w*    2 
V  (inZ) 

.01004 

.00976 

.00999 

.00970 

.00985 

.01016 

.00990 

.00988 

F. 

l 

0.994 

0.965 

0.989 

0.960 

0.975 

1.006 

0.979 

0.978 

Win"4) 

.1128 

.1060 

.1116 
avg. 

.1048 
.1096 

.1082 

.1156 

.1093 

.1090 
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TABLE  IV 
FLOW  DISTRIBUTION  THROUGH  ORIFICE  PLATE, 
LOSS  COEFFICIENT  PLANE  D  TO  PLANE  X 

RUNS  R4,  R5,  R6 

Location 

1 

2 

3 

4 

5 

6 

7 

8 

R4      w  *  =  .61799,  K  =  .79,  P.   =  P_  -  .508  A  P_   =  1122.90 
D                       X 

P  -P 
D  E   ^x 

(Cm  H^) 

104.70 

99.50 

102.40 

96.60 

101.60 

106.70 

101.50 

100.50 

Re(l0-4) 

d=0.3 

2.752 

2.687 

2.724 

2.650 

2.711 

2.780 

2.711 

2.687 

w*    2 
D.  (in  ) 

.01188 

.01160 

.01176 

.01144 

.01170 

.01200 

.01170 

.01160 

F. 

l 

1.000 

0.976 

0.989 

0.962 

0.986 

1.008 

0.985 

0.981 

Vx(i""4) 

.1185 

.1127 

.1160 
avg. 

.1094 
.1151 

.1150 

.1208 

.1149 

.1138 

R5      w  *  =  .68534,  K  =  .79,  P^  =  P  -  .50  A  P.   =  1152.88 

A 

D  E 
(Cm  HP) 

140.95 

135.05 

137.25 

132.05 

137.95 

142.95 

136.45 

135.95 

Re  (ICT4) 
d=0.3 

3.200 

3.139 

3.158 

3.106 

3.168 

3.220 

3.146 

3.146 

*    2 
V  (in  ) 

.01339 

.01313 

.0321 

.01299 

.01325 

.01347 

.01316 

.01316 

F. 

l 

1.016 

0.996 

1.003 

0.986 

1.005 

1.022 

0.999 

0.999 

kD-X(i""4» 

.1228 

.1177 

.1196 

.1151 

.1202 

.1246 

.1189 

.1185 

R6      w   =  .75132,  K  ■  .79,  P   =  P  -  . 50  A  P    =  1181.06 

A 

p  -P 
D(Cm  H20) 

182.90 

176.05 

176.90 

171.35 

180.90 

185.60 

177.90 

177.05 

Re(IO-4) 
d  =  0.3 

3.681 

3.618 

3.626 

3.574 

3.661 

3.706 

3.606 

3.628 

* 
w^   , .  2 . 
Di  (in  ) 

.01479 

.01454 

.01457 

.01436 

.01471 

.01489 

.01461 

.01458 

F. 

l 

1.024 

1.006 

1.008 

0.994 

1.019 

1.030 

1.011 

1.009 

Vx'"1"' 

.1275 

.1228 

.1233 

av£. 

.1195 
.1245 

.1261 

.1294 

.1240 

.1234 
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TABLE  V 
AREA  RESTRICTION  PLANE  E,  MACH  #  PLANES  E,  F  &  I, 
LOSS  COEFFICIENTS  PLANES  D  TO  F ,  X  TO  F ,  F  TO  G  &  G  TO  I 

RUN 

Rl 

R2 

R3 

R4 

R5 

R6 

E   (cm  H20) 

1169.33 

1159.50 

1096.98 

1125.75 

1155.63 

1184.63 

PE    (cm  H20) 

1086.43 

1082.89 

1061.75 

1072.87 

1084.22 

1091.77 

\ 

.32 

.31 

.21 

.26 

.30 

.34 

TD 

.979 

.981 

.991 

.986 

.982 

.977 

—  * 

WE 

ST 

.81367 

.78919 

.55660 

.67181 

.76516 

.86410 

(£  E 

.567 

.571 

.588 

.582 

.574 

.570 

PF    (cm  H20) 

1148.37 

1140.39 

1087.82 

111.81 

1137.18 

1161.11 

Pt 

F    (cm  H20) 

1161.24 

1152.46 

1093.46 

1120.18 

1148.11 

1174.90 

% 

.12 

.12 

.08 

.10 

.12 

.13 

k          -4 
D-F    (in   ) 

.1309 

.1285 

.1166 

.1212 

.1280 

.1331 

kx.F    (in"4) 

.0070 

.0068 

.0015 

.0061 

.0083 

.0086 

VG    (in_4) 

.0293 

.0288 

.0261 

.0284 

.0290 

.0295 

PI     (cm  H20) 

408.34 

407.81 

407.71 

405.50 

Pt 

I    (cm  H20) 

420.92 

426.43 

432.78 

437.57 

*I 

.21 

.25 

.29 

.33 

k           -4 
Vl    (in   ) 

.0492 

.0501 

.0504 

.0517 
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TABLE    3ZI 
PERIPHERAL   DISTRIBUTION   OF   LOSS   COEFFICIENT 
PLANES  A  TO  F    AND    PLANES    D    TO    F 

KD-F 

RUN 

RedO"6) 

F2 

F4 

F6 

F8 

Rl 

1.307 

0.986 

0.989 

1.038 

1.001 

R2 

1.262 

0.987 

0.991 

1.017 

1.004 

R3 

0.868 

0.967 

0.997 
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0.987 

KD-F 

R4 

1.028 

0.978 

1.002 

1.017 

1.002 

R5 

1.177 
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1.003 
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RIO 

1.276 

0.997 
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Rll 
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1.005 
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RESULTS 

TABLE   IX 
OF  WATER  TESTS  WITH  SCALE 
MODEL  NOZZLE   TUBE 

RUN 

DATA 
POINT 

W 
1 b/sec. 

in3/ sec. 

Pin2 
lb/in2 
gage 

Afv  -■ 

'in-dis. 

lb/in2 

RedO*4) 
d=.3l6 

^in-dis. 

WTl 

1 

0.384 

10.639 

9.175 

9.624 

2.67 

11.187 

t 

•H 
CO 

a 

m 
oo 
r^ 

o- 

■—1 

II 

B 
u 
cfl 

PL, 

2 

0.495 

13.714 

14.875 

15.621 

3.44 

10.928 

3 

0.578 

16.014 

20.625 

21.642 

4.03 

11.104 

4 

0.626 

17.344 

24.200 

25.388 

4.36 

11.105 

5 

0.656 

18.169 

26.000 

27.304 

4.56 

10.883 

6 

0.693 

19.194 

27.950 

29.411 

4.82 

10.504 

7 

0.782 

21.681 

34.150 

36.014 

5.45 

10.081 

WT2 

1 

0.254 

7.037 

3.500 

3.696 

1.77 

9.822 

1 

i 

2 

0.238 

6.594 

3.000 

3.172 

1.65 

9.600 

3 

0.382 

10.584 

8.675 

9.119 

2.66 

10.711 

4 

0.448 

12.412 

11.950 

12.560 

3.12 

10.728 

5 

0.521 

14.435 

16.300 

17.126 

3.63 

10.815 

•H 

a 
o 

00 

r- 

<t 

i— i 

II 

E 

4-1 

6 

0.560 

15.515 

18.475 

19.430 

3.90 

10.620 

7 

0.612 

16.956 

22.000 

23.340 

4.26 

10.682 

8 

0.649 

17.981 

24.550 

25.832 

4.51 

10.513 

9 

0.690 

19.117 

27.375 

28.824 

4.80 

10.378 

P 

1 

f 

10 

0.732 

20.281 

30.225 

31.856 

5.10 

10.191 

11 

0.770 

21.334 

32.950 

34.755 

5.35 

10.047 

12 

0.809 

22.414 

35.150 

37.142 

5.64 

9.728 

13 

0.823 

22.802 

37.350 

39.412 

5.74 

9.974 

14 

0.870 

24.104 

39.600 

41.903 

6.06 

9.490 

15 

0.887 

24.575 

41.150 

45.545 

6.17 

9.487 
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TABLE   X 

RESULTS   OF  WATER  TESTS  WITH  GEOMETRICALLY   SIMILAR 

MODEL  OF   SECTION  BETWEEN  PLANES  W  AND  X 

RUN 

W 

lb/ sec 

inVsec 

Re(lO"4) 
d  =  1.05 

U 

*.-. 

K(.757) 
pos.  B 

K(76l) 
pos.  A 

KC7I4) 
pos.  A 

K(.7I4) 
pos.  B 

Wl 

0.198 

5.488 

0.415 

1.082 

5.824 

.720 

.697 

0.313 

8.674 

0.656 

1.388 

6.886 

.768 

.727 

0.509 

14.112 

1.068 

1.167 

5.600 

.775 

.739 

0.739 

20.474 

1.549 

1.253 

5.304 

.787 

.787 

0.910 

25.203 

1.907 

1.232 

5.333 

.780 

.778 

1.007 

27.912 

2.112 

1.219 

5.375 

.774 

.771 

1.124 

31.149 

2.357 

1.219 

5.419 

.780 

.768 

1.191 

32.994 

2.496 

1.204 

5.402 

.766 

.764 

1.270 

35.203 

2.663 

1.203 

5.367 

.767 

.767 

1.372 

38.023 

2.877 

1.204 

5.576 

.755 

.753 

1.515 

41.969 

3.175 

1.208 

5.673 

.743 

.748 

1.759 

48.715 

3.686 

1.216 

5.689 

.745 

.749 

1.829 

50.682 

3.834 

1.222 

5.680 

.746 

.751 

W2 

0.878 

24.340 

1.841 

1.224 

5.664 

.762 

.753 

1.183 

32.769 

2.479 

1.208 

5.617 

.760 

.751 

1.344 

37.246 

2.818 

1.181 

5.718 

.752 

.736 

1.437 

39.814 

3.012 

1.178 

5.730 

.745 

.740 

1.812 

50.213 

3.799 

1.174 

5.619 

.754 

.740 

1.929 

53.435 

4.044 

1.207 

5.686 

.747 

.746 

1.946 

53.913 

4.078 

1.193 

5.768 

.742 

.737 

W3 

0.559 

15.485 

1.221 

1.487 

8.669 

.794 

.748 

0.809 

22.406 

1.767 

1.481 

8.456 

.802 

.748 

0.876 

24.285 

1.915 

1.481 

8.688 

.792 

.734 

0.988 

27.384 

2.160 

1.440 

8.418 

.793 

.746 

1.078 

29.881 

2.357 

1.462 

8.693 

.786 

.730 

1.262 

34.958 

2.757 

1.457 

8.971 

.773 

.724 

1.343 

37.222 

2.936 

1.454 

9.112 

.766 

.715 

1.440 

39.903 

3.147 

1.449 

9.035 

.768 

.718 

1.494 

41.407 

3.266 

1.452 

9.230 

.760 

.714 

1.589 

44.029 

3.473 

1.447 

9.034 

.767 

.724 

1.690 

46.828 

3.693 

1.434 

9.216 

.756 

.713 

1.754 

48.587 

3.832 

1.430 

8.698 

.778 

.733 
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TABLE  XI 
OXIDIZER  DISCHARGE  &  RETURN  ASSEMBLY,  INLET 
ORIFICE  DATA  AND  EQUIVALENT  FLOW  RATES 

P  fc   -  14.743  lb/in2 

RUN 
AVERAGE  OF  DATA  POINTS 

R7 

R8 

R9 

R10 

Rll 

P     (lb/in2) 
u 

18.055 

18.943 

19.986 

21.179 

22.220 

Ah    (in  Ho0) 
o        2 

0.991 

1.198 

1.448 

1.736 

1.951 

T     (°R) 
u 

541.0 

564.5 

561.0 

554.0 

580.5 

Ven  <°R> 

539.5 

559.5 

559.0 

553.5 

576.0 

Re<10"5>  d=4.07" 

0.917 

0.985 

1.117 

1.276 

1.320 

*       2 
w     (in  ) 
u     N    ' 

.50180 

.53848 

.57723 

.61355 

. 63494 

P      (lb/in2) 
Z 

17.769 

18.651 

19.602 

20.750 

21.721 

Pt     (lb/in2) 

A 

17.597 

18.388 

19.296 

20.393 

21.330 

Pt     (lb/in2) 
F 

15.589 

15.827 

16.114 

16.462 

16.752 

Pt     (lb/in2) 
G 

15.464 

15.668 

15.913 

16.216 

16.463 

-*      „  2. 

«z     (m  ) 

.50916 

. 54449 

.58750 

.62596 

.64702 

—  *       2 
WA     (in  ) 

.51414 

.55228 

.59679 

.63690 

. 65888 

wp     (in  ) 

.58037 

.64164 

.71467 

. 78901 

.83895 

2 
«G     (in  ) 

.58505 

.64817 

.72367 

.80098 

.85364 

?t     (lb/in2) 
o 

18.005 

18.890 

19.894 

21.093 

22.132 
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TABLE  XII 
LOSS  COEFFICIENTS  PLANES  Z  TO  A ,  A  TO  F , 
F  TO  G,  G  TO  H,  H  TO  I ,  W  TO  X 
DISCHARGE  MACH  NUMBER 

RUN 

R7 

R8 

R9 

RIO 

Rll 

Va     <in_4> 

.0374 

.0473 

.0451 

.0438 

.0430 

Z                               f      "4N 

kA-F      (Xn   > 

.4316 

.4564 

.4630 

.4752 

.4944 

Vg      (in_4) 

.0237 

.0245 

.0243 

.0240 

.0244 

kG-H      (in"4) 

.0414 

.0415 

.0410 

.0407 

.0409 

P       (in  HO) 
I 

420.98 

425.07 

429.52 

434.54 

439.60 

Px      (in  H20) 

407.00 

407.03 

406.71 

406.23 

405.83 

\ 

.22 

.25 

.28 

.31 

.34 

kG-I      (ln   > 

.0500 

.0486 

.0486 

.0506 

.0493 

Vi        (in"4) 

.0086 

.0071 

.0076 

.0099 

.0084 

*A-X      (in_4) 

.4296 

.4539 

.4600 

.4702 

.4894 

P       (lb/in2) 
X 

15.601 

15.845 

16.138 

16.492 

16.801 

Re  <10'4)  d  =  0.3 

2.392 

2.567 

2.913 

3.329 

3.442 

P       (lb/in2) 

17.338 

18.024 

18.914 

19.935 

20.752 

— *          2 
«w        (^  ) 

.52182 

.56342 

.60886 

.65180 

.67614 

J  =Apt      /Apb.c 

'w-x/ 

1.46 

1.458 

1.453 

1.447 

1.445 

Vx    (in_4) 

.3678 

.3908 

.3959 

.4061 

.4193 

£w-x 

9.938 

10.560 

10.698 

10.973 

11.329 
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APPENDIX  A 
DERIVATION  OF  FLOW  MEASURING  EQUATION 

The  pressure  drop  P  -  P,  across  the  flow  nozzle  is  referred  to 

the  absolute  pressure  P  upstream  and  is  to  enter  calculations  in  the 
r         u 

Pu-Pd 
dimensionless  form  — .   Thus  any  consistent  unit  of  pressure  may  be 


used. 


The  equation  for  flow  rate  w  (lb  /sec)   is 


=    •*V,A.f4l  fr   ftffr  (V  Pd) 


where: 

*£■     -  nozzle  discharge  coefficient  (independent  of  Re) 

Y,  -  expansion  factor  to  account  for  compressibility  effects 

2  2 

A,  -  area  based  on  nozzle  diameter  (in.  )  ■  7.17263  in. 

f  -  area  multiplier  to  account  for  the  thermal  expansion  of  the 

nozzle 

f   -  Reynolds  number  correction  factor 

g  -  32.17  ft/sec2 

o 

Pu  -  upstream  density  (lb  /ft  ) 

m 

2 
Pu  -  pressure  upstream  (lb/in.  ) 

2 
Pd  -  pressure  downstream  (lb/in.  ) 

Using  the  equation  of  state  for  a  perfect  gas 


V  u     


(A2) 


*TU 


VDI  Durchfluss  Messregeln.  DIN  1952,  VDI  Verlag,  Berlin, 
1943,  p.  5. 
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where  R  ■  53.35  ft-lb/lb  R,  and  T  is  the  upstream  temperature  (  R) , 

m        u  N 

equation  (Al)  can  be  rearranged  to  give  the  equivalent  flow  rate 


^±  J/!  -       y,  Ad  ft  f'     2  )/Pu  -  P„    M)     (A3) 


The  dimensionless  coefficients  must  be  investigated  for  the  noz- 

zle  with  diameter  ratio  6  *  0.7434  and  m  =  ft   =  0.552  where  6  =  d/D, 

2 
The  discharge  coefficient   oC  for  m  ■  0.552  is  found  to  be  1.114. 

3 
The  density  coefficient  Y,  is  obtained  from  the  ASME   standards. 


Fo 


Pu-Pd 


r  air  with  <f  =  c  /c  =  1.4  and  B  ■  .743  using  x  ■  — r there  is 

p     v  V  6  Pu 

/  /.+l&  0.2.86 

/  X  /    /   -    0.30S     (l-K)     +         J 


4 
For   stainless   steel   the  area  multiplier   f     can  be  expressed  by 


f       =         1.0        +      0.00  IS     f    ^u  '    6Q  \ 
i  l        /o©        1 


(A5) 


The  design  of  the  nozzle  is  such  that  f  =  1  for  all  pressure  ratios. 
Combining  the  above  quantities  equation  (A3)  becomes 


2 
Ibid.  Appendix  6. 

3 
Supplement  to  ASME  Power  Test  Codes,  Chap.  4,  Part  5,  ASME, 

New  York,  1959,  p.  76. 

4 
Stearnes,  R.  F.  et.  al.  Flow  Measurements  with  Orifice  Meters, 

Van  Nostrand  Co.,  Inc.,  New  York,  1951,  p.  257. 
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APPENDIX  B 
DATA  REDUCTION  FOR  AIR  AND  WATER  TESTS  WITH  FLOW  CHANNEL 

B.  1  Air  Tests 

A  total  of  seven  data  points  was  observed.   The  quantities  mea- 
sured were: 

H  -  pressure  upstream  of  flow  measuring  orifice  against  atmo- 
sphere (in-Hg) 


Ah 


-  change  in  static  pressure  across  orifice  (in.H«0) 

t   -  temperature  upstream  of  orifice  (  F) 

t  n    -  temperature  in  plenum  2  (  F) 
plen 

H    -  manometer  reference  (cm  H„0) 
atm  2 

H.   -  inlet  total  pressure  at  Kiel  probe  against  atmosphere 


(cm  H20) 


H , .   -  discharge  total  pressure  against  atmosphere  (cm  H90) 

h    -  30.01  in.Hg 
atm 

The  equivalent  flow  rate  w  is  given  by  equation  (A6  ).   It  is 

convenient  to  use  the  model  inlet  as  a  reference  condition.   Then  the 

equivalent  flow  rate  is 


*        *■       P      /  T 
w  =    w        r*    J     'tin  (in1)       (Bl) 

R    I  T 


where  P  =  P   =  0.49112 
u    t 
u 


{  Hu  +  30.01 i     (lb/in   ) 

P.      =  Pt        =  0.49112    (  0.02894   (H   ,      -   H,    )  +  30.01  I     (lb/in2) 
in  t.  x   atm  in  ) 

in  v 

Reynolds   number  Re   based  on   inlet   diameter  d   is 

Ke    -  w-    =    — — (B2) 
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m 


Equation  (B2)  can  be  rewritten  as 

Re  -  <  -^  J*'1*)—  (.3) 

Air  temperature  for  these  tests  varied  between  85  and  95  F.   Introducing 

-7         2 
a  representative  viscosity  of  3.9  x  10   lb-sec/ft  ,  in  addition  to 

other  non-varying  quantities,  into  equation  (B3)  gives 

Re  --   wi*        _^_    C4.7£e>*/o5)  (B4) 

The  loss  coefficient  is  determined  by 
P      -     P  / 

*,„_*  =  Jl ^/(w.*)  ('*'')  (B5) 

where  P     =  0.49112  \  0.02894  (H    -  H    )  +  30.0l\   (lb/in2) 
t,.  (        v  atm    t, .  ) 

dis  v  dis 

B.2  Water  Tests 

A  total  of  four  runs  was  made,  each  consisting  of  several  data 
points.   The  runs  were  numbered  RMW  1,  etc.   The  quantities  measured 

were: 

2 
P.    -  inlet  static  pressure  (lb/in  gage) 

W  -  weight  of  water  (lb) 

t  -  elapsed  time  (sec) 

Volume  flow  rate  \f   is  obtained  by 

=     ^   POO)  (^/SSO  (B6) 

3.  60  9"5 


Re  based  on  inlet  diameter  is 

/(J,      %  <*  *  (B7) 
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where  V  ■  kinematic  viscosity,  for  water  at  72  F  "   152.257  in  /sec 

Equation  (B7)  is  reduced  to 

Rz    '         ^  /   2-39.  /64-5  (B8) 

Because  the  inlet  and  discharge  areas  are  equal 

APt       -  Ap.   ..   -  P. 
t,   ,.      m-dis    in  g 
in-dis 

The  loss  coefficient  is  determined  by 

*iVdis      * <'*  )   <B9) 

where  V   =  specific  weight  of  water  at  72  F. 

1j 
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APPENDIX  C 
DATA  REDUCTION  FOR  TESTS  WITH  ARES  MODEL 

C.l  General 

Methods  of  data  reduction  used  to  obtain  loss  coefficients  and 
other  meaningful  quantities  from  tests  on  the  ARES  model  (as  described 
in  Section  5  of  the  main  body),  are  discussed  in  the  following  Sections. 
They  are  listed  in  the  same  sequence  as  used  for  the  performance  of  the 
tests. 
C.2  Oxidizer  Return  Sub-assembly 

The  quantities  necessary  to  obtain  mass  flow  rate  are  identical  to 

those  described  in  Appendix  B.   The  equivalent  flow  rate  w  is  given  by 

equation  (A6) .   In  addition  the  following  quantities  were  recorded 

H  ,    -  pressure  in  plenum  2  (cm  H„0) 
plen  2 

H    ,  -  atmospheric  reference,  manometer  bank  1  (cm  H_0) 
atm  1  2 

H.   -  pressure  at  Kiel  probe  (cm  H~0) 

H_  -  pressure  in  plane  D  (cm  H?0) 

h   -  static  pressure  at  selected  locations  in  annulus  (in.Ho0) 
an  x    2 


H 


an 


total  pressure  at  selected  locations  in  annulus  (in.H90) 


H  .   0  -  atmospheric  reference,  manometer  bank  3  (in.Ho0) 
atm  3        r  v    2  ' 

h  -  static  pressure  in  plane  I  (in.H~0) 
H  -  total  pressure  in  plane  I  (in.H^O) 
t  1    -  temperature  in  plenum  2  (  F) 
From  the  polaroid  pictures  the  following  quantities  were  read 

h     -  static  pressure  in  plane  E  at  holes  7,  13,  20,  26,  33,  39, 

46,  52  (cm  H20) 
h  -  static  pressure  in  plane  F  at  holes  1,  9,  17,  25  (cm  Ho0) 
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H  -  total  pressure  in  plane  G  between  holes  32-1,  4-5,  9-10, 
G 

13-14,  17-18,  20-21,  25-26,  29-3  (cm  H20) 

h     -   static  pressure   in  plane  H  above  discharge  holes    11,    14,    18, 
H 

22,   26,    30,    31,    34,    38,   42,   46,   50   (cm  H20) 


Re=     wu    J\_        4  (\t.) w*      P^    o.oiohz     (cl) 


The  loss  coefficient  k  for  flow  between  each  plane  is  referred  to 
an  inlet  condition.   The  equivalent  flow  rate  must  then  be  referred  to 
conditions  at  intermediate  planes.   For  example^an  average  equivalent 
flow  rate  at  plane  F  is 


£  *  *     w  *   ?*        I T0  (C2) 

It  «*        


where  ?fc   -  0.49112  ^0.02894  (H^  \   -   H^l  +  P^  (lb/in2) 
F  # 

T  -  temperature  in  plane  D  assumed  to  be  the  same  as  that  in 

plenum  2 

P  and  T   -  determined  as  in  equation  (Bl). 
u      u 

—  *  —  * 
In  a  similar  fashion  the  equivalent  flow  rates  w  ,  w   are  determined. 

A  geometrically  similar  model  of  a  tube  section  from  plane  D  to 
plane  X  was  tested  by  Vavra.   Data  from  these  tests  plotted  in  Figs.  29 
and  30  include  the  orifice  discharge  coefficient  K,  and  the  loss  in 
total  pressure  "f     from  plane  D  to  plane  X  divided  by  the  static 
pressure  change  across  the  orifice.   For  each  orifice  preceding  a  mea- 
suring tap  in  plane  E,  an  equivalent  flow  rate  is  found  by 
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# 


«o.     -     *Y.A.      z     IP.-  p.  (c3) 


l 


2 
where  A,  =0.0363  in  (area  of  orifice  opening). 

P„  and  P^  are  determined  as  was  P^  above.   For  air  with 
ED  tp 

>  ■  C  /c  *  1.4  from 
P  v 

P  -  P 
V   a    /.0   -  0.3  6     *>    E  (C4) 


Equation  (C3)  becomes 


* 


P.-  P. 


W    =   O.OS204-  /<  Y.  i  /  rp "   £ 
O.  '  1/  —^ <C5) 


L 


PD 


Equation  (B4)  is  modified  to  find  Re  referred  to  tube  inlet  diameter 


*    P  7 


Re  ~  "D.    — ^_   V  ■"***  ^  »«  '  (C6) 

This  relation  contains  a  representative  value  of  viscosity  for  the  tem- 
perature range  encountered  in  the  tests.   An  iteration  is  carried  out 

using  equations  (C5)  and  (C6) .   A  value  of  K  is  assumed  and  equation 

* 
(C5)  is  solved.   Using  this  value  of  w   ,  Re  is  calculated  from  equa- 

i 

tion  (C6) .   A  corrected  value  of  K  is  then  read  from  Fig.  29  and  the 

process  continued.   Generally,  one  iteration  is  required  only. 

The  peripheral  flow  distribution  at  the  orifice  plate  is  found 

from 

F,  «   (£2.)      wo,  /-  *  (C7) 
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Tht  total  pressure  P   is  found  from 

X 


where  \    is  read  from  Fig.  30  for  the  Re  found  with  equation  (C6) .   It 
is  possible  also  to  estimate  P   by  using  equation  (C7)  with  t  =0.48, 
a  standard  value  for  sharp  edge  orifices  with  p  =  0.72. 

The  flow  area  in  plane  E  is  3.676  sq.  in.   It  is  possible  to  find 
an  average  area  restriction  factor  $„.   The  equation  of  continuity  for 
average  conditions  is 

V  =■    (C9) 

I    A    0 

Also 

-  Sl       =    5  (  P.  -  P  \ 


Combining  equations  (C9)  and  (CIO)  gives 


Vs/ 


i/»a#x 


Zf[~P<-P\  (cu) 


Using  the  equation  of  state  for  a  perfect  gas  P  ■  p  RT,  equation  (Cll) 
becomes 


P  2  J  1*1       *      $ 


(C12) 
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Referring  equation  (C12)  to  static  conditions  at  plane  E  and  solving 
for   <£ 

v2 


*£  = 


L     ("fci  ) 


P 


(C13) 


where 


*•        *   p 


Vt. 


An  average  Mach  number  M  is  determined  from  the  pressure  ratio 

P _/P  .   Assuming  no  change  in  total  temperature  across  the  orifice, 
E   t_ 
E 

i.e. ;  T   ■  T  ,  it  is  possible  to  find  the  temperature  ratio  needed  to 

E 

solve  equation  (C13). 

With  the  equivalent  flow  rates  and  equation  (12)  ,  average  loss 

coefficients  1c   ,  k   ,  k ,  k    and  k    are  determined.   For  example 
D~F    F~G    G~H   G~l      D~X 

the  average  loss  coefficient  from  plane  F  to  plane  G  is 


k 


Pt.-  p* 


F-fi 


<5  2 


a«~4) 


(C14) 


At  a  particular  location  in  plane  I  the  discharge  velocity  is  pro- 
portional to  the  square  root  of  the  dynamic  pressure. 


i 

V,  = 


H*t„3-  M;]'Kw  M=   ?:    <''*•"*'  (ci5) 


At  each  location  the  velocity  relative  to  the  average  is  then 
V:/_   _ 


V 


V  *  ?  v^ 


(C16) 
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C. 3  Scale  Model  Nozzle  Tubes 

The  following  quantities  were  recorded 
t'  -  elapsed  time  (sec) 

W  -  weight  of  water  collected  (lb) 

2 
P^    -  inlet  static  pressure  (lb/in  gage) 
in  g 

t  -  water  temperature  (  F) 

2 
A,   -  inlet  area  (0.07843  in  ) 
in 

2 
A,.   -  discharge  area  (0.11341  in  ) 

h  ,      -   barometric  pressure  (in  Hg) 
atm  r        \  or 

Flow  rate  w  was  determined  by  dividing  the  weight  of  water  by  the 
elapsed  time.   More  useful  for  calculations  is  the  volume  flow  rate 


3.6095 


(C17) 


Velocities  at  inlet  and  discharge  can  be  written  as 

V  =   V>/4  (C18) 

Assuming  atmospheric  pressure  at  the  discharge  the  loss  in  total 
pressure  through  the  tube  is 


p^ '  p'dl5  "  i p-  "  P*^ +  4  f4 (v-  "  < ]      <  "■  '*> 


(C19) 


In  all   tests    the  water   temperature  varied  only   from  55      to   59   F.      For 
t  ■  57   F  equation    (C19)   becomes 

****.««,      =      P>°>       ¥       (3-9&S    *     ">'''*'  (C20) 

Re  referred  to  an  inlet  diameter  of  0.316  inches  is  obtained  from 
equation  (B7)  with  V   -  160.279  x  10"5  in2/sec 

/?e  -    ^X/Q*  (C21) 

3.9  8 
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Using  the  inlet  as  reference,  the  loss  coefficient  F  is 


in-  c<i«» 


(C22) 


V* 


-   0.13S   x  /o'4       /6-sec* 


y 


i*  + 


C.4  Geometrically  Similar  Model  of  Section  from  Plane  W  to  Plane  X 

Figure  28  shows  6  manometer  tubes  in  a  schematic  manner.   Actually, 

Stations  2  through  6  have  2  taps  each  so  that  h»  (in-H-O)  is  an  average 

pressure  at  Station  2.   In  addition  to  pressures,  the  mass  flow  rate 

and  water  temperature  were  recorded  as  described  in  sub-section  C.3. 

Volume  flow  rate  is  determined  as  in  the  former  tests.   Reynolds  number 

referred  to  an  inlet  diameter  of  1.05  in.  is 

p       \f  x  /o  + 
13.  21  S 

for  average  temperatue  conditions. 

It  is  assumed  that  taps  3  and  5  are  placed  such  that  the  flow  has 
the  same  area  as  the  orifice  opening.  The  error  in  this  assumption  is 
absorbed  in  the  flow  coefficient  K. 

Across  an  orifice  the  equation  of  continuity  gives 

^d  ^d  =  Kr  Vor  (C24) 

where  d  -  station  ahead  of  orifice 

or  -  station  behind  the  orifice 
Expressing  terms  in  in.  of  liquid  head,  Bernouilli's  equation,  ignoring 
losses,  gives 
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23  2* 


Combining  equations  (C24)  and  (C25) 


v._  = 


or 


2  q,  (  Kj  -  Kor) 


i 


I- 


r 


(C26) 


If  losses  and  other  factors  are  included  in  a  coefficient  K,  the  volume 
flow  rate  can  be  written  as 

V  -    */<       ■     K  A.r  Vor  (c27) 

Equations  (C26)  and  (C27)  are  combined  to  give 


(C28) 


«     d-or 

Three  orifice  sizes  were  tested.   Table  C-I  lists  factors  for  each 
orifice. 

TABLE  C-I 
GEOMETRIC  FACTORS  FOR  VARIOUS  ORIFICES 

DIAMETER        6  1-  fr>         AQr 

0.757      0.7209      0.7298      0.45007 
0.761      0.7248      0.7241      0.45484 
0.714       0.6800      0.7862      0.40039 
Using  Station  6  as  an  example  and  assuming  stagnation  conditions 
at  Station  1,  the  head  loss  IL  is  written  as 

M,-"4  -   V*4-   &'  -        "l/2  (C29) 
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Compared  with  other  losses,  the  loss  through  the  inlet  orifice  is 
negligible.   A  loss  coefficient  from  the  beginning  of  the  1.05  in.  pipe 
to  Station  6  is  then 


(C30) 


?,t 

= 

A 

f 

*4 

Introducing 

the 

volume   flow 

rate 

6 

- 

578 

,  896 

(C31) 

C.5  Discharge  and  Return  Assembly 

In  addition  to  the  quantities  measured  in  the  return  sub-assembly 

the  following  pressures  were  read 

h  -  static  pressure  at  entrance  of  model  (cm  H~0) 

H  -  total  pressure  at  entrance  of  model  (cm  H«0) 

h    „  -  atmospheric  reference,  manometer  bank  2  (cm  Ho0) 
atm  L  i 

h  -  static  pressure  in  plane  A,  6  equally  spaced  taps  (cm  H„0) 

IL  -  total  pressure  in  plane  Z  (cm  H?0) 

h  -  static  pressure  in  plane  A  at  holes  1,  14,  27,  40   (cm  H90) 

H  -  total  pressure  in  plane  A  between  holes  7-8,  20-21,  33-34, 

46-47  (cm  H20) 

h  -  static  pressure  in  plane  B  at  holes  7,  13,  20,  26,  33,  39, 
B 

46,  52  (cm  H20) 
h  -  static  pressure  in  plane  C  at  holes  7,  13,  20,  26,  33,  39, 

46,  52  (cm  H20) 
H  1    -  pressure  in  plenum  2  referenced  to  atmosphere (cm  Hg) 
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■ff      -k 

Equivalent  flow  rates  found  as  in  Section  C.2  include  w  ,  w  , 
^  u    z 

— —  "M      ~~  *Jc  ***   *te 

w   ,  w   ,  and  w  .   These  flow  rates  are  used  to  calculate  the  loss 
A  *   F       G 

coefficients  \_^   k"A_p,  kp_G>  k^,  and  k^. 

Using  Re  with  the  4  in.  diameter  of  the  supply  pipe  as  character- 
istic length,  the  argument,  k^  is  read  from  Fig.  35.  Total  pressure 
at  plane  X  is  then  given  by 

where  :  —  — 

KAOC  Vf         kx-f 


K      -        ^       -     f   kA  h      W  *  Z    I  (C32) 


4-X 
An  average  Re  of  the  flow  from  plane  W  to  plane  X  can  be  written  as 


Re  -  ™  v      t-w-r)  (c33) 


"W     FF 


where: 

* 
W 

K 

l/T» 

pt 

V  r 

xw 

1/               IA. 

(C34) 


and 


\,  =   f  ^  P.  J  *   Ptx  (C35) 


t 


An  iteration  of  equations  (C33) ,  (C34) ,  and  (C35)  is  carried  out  by 
assuming  an  initial  value  of  T  from  Fig.  46.   The  process  converges 

rapidly  because  £  has  a  limited  range  of  values. 

—      —  *  — 

With  the  computed  values  of  P   and  w„  ,  the  coefficient  k_  _  is 

S,      W  "D-X 

calculated.   From  equation  (12) 

^W-X  ■    *W'A  (C36) 

24* 
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